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ABSTRACT 


Undisturbed  coal— bear! ng  strata  of  the  Interior  Plains 
region  are  commonly  covered  to  a  varying  depth  by  pre- 
glacial  and  glacial  deposits  which  obscure  the  coal  zones 
from  direct  observation*  Indirect  methods  such  as  use  of 
exploration  boreholes  result  in  a  fragmented  view  of  the 
coal  deposit  and  hence  uncertainty  in  the  estimate  of 
distribution  and  quantity  of  coal  resource. 


Statistical  methods  provide  quantitative  statements  of 
geologic  assurance  for  estimators.  Kriging  and  the 
arithmetic  average  are  unbiased  estimators  of  the  coal  zone 
average  thickness  over  areas  sampled  Irregularly  and 
systematically  respectively.  Subjective  probability  methods 
and  an  approximation  formula  provide  uncertainty  estimates 
for  the  area  underlain  by  coal. 

Application  of  these  methods  to  data  from  three  coal 
zones  in  the  Tertiary  Ravenscrag  Formation  of  southern 
Saskatchewan  results  in  estimates  classified  as  measured  and 
indicated  according  to  traditional  resource  classification 
schemes.  In  all  cases  examined*  over  two— thirds  of  the 
uncertainty  in  estimates  of  volume  is  contributed  by  the 
area  estimate  whereas  uncertainty  in  the  average  thickness 
of  the  coal  zone  accounts  for  less  than  one— third. 
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glacial  action  or  salt  collapse  origin  impairs  prediction  of 
the  coal  subcrop  location.  The  irregular  nature  of 
lithologic  boundaries  between  coal  and  penecont emporaneous 
clastic  facies  forms  the  other  major  component  of 
uncertainty  in  the  area. 

The  precision  of  average  thickness  estimates  is 
controlled  by  the  variability  of  the  selected  net  coal 
thickness*  the  size  of  the  area  estimated  and  the  average 
spacing  of  observations.  For  this  reason  classification 
methods  based  only  on  arbitrary  minimum  observation  spacing 
fail  to  classify  accurately  resource  estimates  of  different 
coal  zones  and  different  sizes  of  areas.  A  systematic 
sampling  scheme  on  a  square  grid  is  recommended  as  the  most 
efficient  form  of  exploration  for  resource  evaluation  of 
undisturbed  coal— bearing  strata  in  southern  Saskatchewan. 
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CHAPTER  1 


1  •  Introduct ion 


Decision  makers  in  industry  an 
estimates  of  the  quantities  of  coai  res 
region  or  country  as  basic  information 
among  various  energy  alternatives 
estimates  should  include  some  quanti 
assurance  that  the  resource  is  actually 


The  quantity  of  coal  resource  co 
is  simply  the  product  of  the  area  under 
the  average  thickness  of  the  bed  and  a 
the  boundaries  of  such  a  coal  bed  we 


and  the  bed  was  of  uniform  thickness 


single  areay  thickness  and  quality 
suffice  for  a  resource  quantity  es 
howevery  there  is  significant  uncerta 


values  of  the  areay  thickness  and  quali 


The  area  underlain  by  a  coal  bed 
certainty  because  erosional  boundari 
transitions  are  commonly  obscured  by  mo 
Using  lithologic  logs  and  geophysic 
boreholesy  the  geologist  creates  struct 
seam  surfaces  and  of  the  erosional  unc 
the  seams*  Using  these  structural 


d  government  require 
ources  available  in  a 
for  making  choices 
•  Fur thermorey  such 
tative  statement  of 
available* 

ntained  in  a  coal  bed 
lain  by  the  coal  bedy 
tonnage  factor*  If 
re  exposed  in  outcrop 
and  qualityf  then 
measurements  would 
tiasate*  In  practice 
inty  about  the  actual 
ty  of  coal  beds* 

is  not  known  with 
es  and  lateral  facies 
re  recent  deposits* 
al  logs  from  rotary 
ural  models  of  the 
onforraities  affecting 
interpre ta tions y  the 
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subcrop  boundaries  of  the  seam  can  be  estimated*  Commonly 
without  the  hard  sediment ologlc a 1  information  found  in 
outcrop  or  continuous  cores*  the  geologist  interprets  the 
available  data  to  produce  a  generalized  facies  model 
outlining  areas  of  coal  non-deposition. 

The  thickness  of  coal  in  a  succession  of  beds  varies 
locally*  Sock  partings  and  high-ash  coal  lenses  which 
separate  beds  of  clean  combustible  coal  result  from  local 
depositional  conditions  which  are  difficult  to  predict  given 
the  usual  scale  of  sample  aval lab li li ty  in  a  regional 
resource  evaluation  program*  Furt hermor e *  control  points  are 
irregularly  spaced  and  clustered  in  locations  of  complex 
geology*  resulting  in  a  sample  set  which  may  not  be 
representative  of  the  whole  area* 

Coal  quality  parameters  such  as  specific  gravity, 
heating  value*  ash*  volatile  and  sulfur  content  must  be 
measured  and  used  in  the  selection  of  local  potentially 
mineable  areas*  Specific  gravity  is  used  in  calculating  a 
tonnage  factor  for  the  regional  resource  evaluation* 


Uncertainty  in  the  geological  interpretation  of 
exploration  data  is  compounded  by  errors  in  measurement  of 
coal  thickness  from  geophysical  logs*  errors  in  determining 
location  ard  elevation  of  boreholes  or  other  sample 
locations  from  topographic  maps  and  analytical  errors  in  the 
determination  of  specific  gravity  and  other  quality  data* 
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How  is  the  exploration  geologist  to  provide  an  estimate  with 
some  quantitative  statement  of  geological  assurance? 


Estimation  methods  have  been  developed  which  permit  one 
to  make  statistical  inferences  about  the  values  of 
population  parameters*  Some  methodSf  although  intuitively 
simple*  are  mathematically  complex  and  require  large  numbers 
of  calculations;  until  the  advent  of  high  speed  computers 
the  cost  of  using  such  methods  was  prohibitive*  Advanced 
estimation  methods  in  combination  with  a  computer 
exploration  data  base  can  nov  provide  ccal  resource  quantity 
estimates  more  rapidly  and  economically  than  previous  manual 
methods* 


This  thesis  deals  with  the  application  of  several 
statistical  methods  to  the  problem  of  ccal  resource  quantity 
estimation  in  the  undisturbed  coal— bearing  Tertiary  strata 
of  southern  Saskatchewan* 


1*1  Bgsgiiffi.fi  Estima.tign  Igrmlpg.lQgy 

* 

McKelvey  (1976)  has  classified  the  resource  base  of  a 
mineral  commodify  along  two  axes  (Figure  1)*  He  restricts 


the  use  of  reserve  to  (p*2): 

"•••that  portion  of  the  identified  resource 
from  which  a  usable  mineral  or  energy 
commodity  can  be  economically  and  legally 
extracted  at  the  time  of  determination*" 

This  study  is  not  concerned  with  the  economics  of  coal 

extraction*  but  rather  it  focuses  on  the  other  axis*  the 
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"degree  of  geologic  assurance".  Therefore  the  quantities 
estimated  are  termed  resources  rather  than  reserves. 

The  degree  of  geologic  assurance  or  "assurance  of 
existence"  is  composed  of  two  separate  problems.  (1)  does 
the  resource  exist?  and  (2)  if  it  exists*  what  is  its  size? 
The  former  question  is  of  paramount  importance  in  estimating 
the  volume  of  undiscovered  oil  and  gas  resources  (Haun, 
1976*  EMRf  1977)*  whereas  the  presence  of  coal  in  the 
Ravenscrag  Formation  of  southern  Saskatchewan  is  well  known 
and  only  the  size  of  the  resource  is  in  question. 

The  degree  of  geologic  assurance  of  coal  resources  has 
been  reported  by  government  agencies  according  to 
classification  schemes  based  on  the  maximum  distance  between 
observation  points.  Table  1  summarizes  these  distance 
cutoffs  in  metres  for  the  "Measured",  "Indicated"  and 
"Inferred"  resourc es/ reserves  classes  of  Canada,  USA,  India 
and  the  USSR.  The  distances  in  Table  1  apply  to  horizontal 
or  gently  dipping  beds  under  conditions  similar  to  those 


found  in  the  Interior  Plains  of  western  Canada. 
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RESOURCES 

jjaasargd 

laf  g.r.r«d 

Canada 

<800 

800-1600 

>1600 

( Plains ) 

USA 

<800 

800-2400 

2400-9600 

India1 

<200 

200-2000 

2000 

USSR1 

365 

365-730 

(uniform  beds) 

275 

275-550 

(fairly  consistent  beds) 

200 

200-365 

(inconsistent  beds) 

1  Reserves 


Table  1:  Coal  resou rce/ r eserve  classifications  according  to 

the  maximum  distance  in  metres  between 
observat ions* 
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The  following  commonly  accepted 
"measured" ,  "indicated"  and  "inferred" 
McKelvey  (  1976*  p*3)* 


definitions  of 
are  taken  from 


it 


Measured 


Resources  are  computed  from  dimensions 
revealed  in  outcrops*  trenches*  mine 
workings  and  drill  holes*  The  points  of 
observation  and  measurement  are  so 
closely  spaced  and  the  thickness  and 
extent  of  coals  so  well  defined  that  the 
tonnage  is  Judged  accurate  within  20 
percent  of  true  tonnage* 


Indicated  —  Resources  are  computed  partly  from 
specified  measurements  and  partly  from 
projection  of  visible  data  for  a 

reasonable  distance  on  the  basis  of 
geologic  evidence* 


Inferred  —  Quantitative  estimates  arc  based  largely 
on  broad  knowledge  of  the  geologic 
character  of  the  bed  or  region  and  where 
few  measurements  of  bed  thickness  are 
available* " 


The  measured  class  is  the  only  one  for  which  the 
uncertainty  is  defined  quantitatively*  The  meanings  of 
inferred  and  Indicated  are  ambiguous  and  hence  subject  to 
misinterpretation*  The  considerable  range  in  distance 
cutoffs  shown  in  Table  1  reflects  varied  definitions  and 
experience  in  resource  estimation*  The  ranges  are  often 
modified  subjectively  to  reflect  local  geology  of  a  coal 
bed*  The  result  of  using  such  classification  schemes  is  an 
inconsistent  partitioning  of  coal  resources  into  classes 
which  have  little  quantitative  value  in  the  decision-making 
process • 


Estimation 


methods 


wh  ich 


pr ovi de 


quan  ti ta t ive 


. 


. 
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statements  of  geologic  assurance  avoid  the  problems  inherent 
in  the  classification  schemes  and  offer  the  advantages  of 
consistency  and  roproduc  ib  il  i  ty  •  Statistical,  estimation 
methods  are  not  necessarily  more  objective  because 
subjective  input  is  required  in  the  geological 
interpretation  of  data,  the  choice  of  a  model  and  the  choice 
of  confidence  level*  Where  data  are  Insufficient  to  permit 
geological  interpretation^  an  estimate  cannot  be  prepared  by 
any  conventional  method  or  statistical  method  based  on 
sampling*  At  low  information  levels,  only  the  methods  of 
subjective  probability  provide  quantitative  statements  of 
assurance • 


Popoff  (1966),  in  a  review  of  conventional  reserve 
estimation  methods,  suggested  a  two— fold  purpose  for  these 
procedures:  ( 1)  valuation  of  the  deposit  and  (21 
illustration  of  the  distribution  of  the  variables*  This 
division  corresponds  to  global  estimation  and  local 
estimation  as  used  by  Matheron  (1971)*  Local  estimates  in 
the  form  of  isopach,  isoquality,  structure  contour  and  other 
isoline  maps  are  used  for  the  selection  of  areas  of  economic 
interest*  Selection  of  coal  which  is  to  be  strip  mined  is 
often  first  performed  using  a  minimum  bed  thickness  and  a 
maximum  overburden  ratio  or  depth*  As  exploration  progresses 
and  information  becomes  available  on  the  hydrogeology  oi  the 
enclosing  beds,  quality  characteristics  of  the  coal  and 


engineering  properties  of  the  enclosing  materials,  further 


. 

. 

- 
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selection  is  made  of  economically  interesting  coal* 

Although  statistical  methods  can  he  and  have  been  used 
for  local  estimation*  the  lack  of  coal  quality  data  and 
engineering  parameters  precludes  the  illustration  of  and 
selection  according  to  these  variables  in  this  study*  Local 
estimation  methods  will  be  used  to  illustrate  the 
distribution  of  variables  and  estimate  subcrop  boundaries 
but  will  not  be  used  for  selection*  Local  estimation  methods 
will  be  discussed  in  the  chapter  on  statistical  methods* 


Resource  estimation  as  used  in  this  study  refers  to  the 
global  estimation  of  the  quantity  of  coal  in— situ  in  a  bed 
without  selecting  local  areas  of  economic  interest* 
Specifically*  the  resource  estimation  equation  is: 

Volume  =  ^rea  *  Jhickness 


1*2  ixsa 

The  objectives  of  this  project  were:  1  )  to  identify 
causes  of  uncertainty  in  the  estimation  of  area  and  average 
thickness  of  coal  zones  using  exploration  data  for  three 
coal  zones  in  two  areas  of  southern  Saskatchewan*  2)  to 
compare  statistical  models  of  the  spatial  variability  of 
seam  thickness  and  structural  configuration  with  conceptual 
models  of  coal  geology*  applying  the  Theory  of  Regionalized 
Variables  developed  by  Matheron  ( 1963)  and  in  particular  to 


apply  the  techniques  of  structural  analysis  and  kriging  to 


. 


' 


the  comparison* 


3)  to  investigate 


alternative  methods  of 


statistics  and  probability  which  can  be  employed  in  coal 
resource  estimation*  and  4)  to  apply  some  of  the  statistical 
methods  to  exploration  data  from  southern  Saskatchewan* 


. 


CHAPTER  2 


2.  gaalflgy  ax  Sglsc.tgri  Sa-al  Zanss 

The  Savenscrag  Formation  of  Paleocene  age  covers  about 
10,000  square  miles  of  southern  Saskatchewan)  and  is  coal¬ 
bearing  throughout  this  area*  Exploration  data  for  three 
coal  zones  in  two  localities  were  selected;  two  in  the 
Estevan  coalfield  (shown  in  Figure  2)  and  a  third  not 
located  in  the  Estevan  coalfield,  chosen  because  of  the 
closely-spaced,  systematic  drilling  pattern*  Due  to  the 
confidential  nature  of  the  data,  the  precise  location  of  the 
third  site  cannot  be  given* 


2* 1  Exgylgas  Igrfr 


Early  work  on  the  geology  of  the  coal— bearing  strata  of 
southern  Saskatchewan  was  restricted  to  exposures  in  the 
Souris  River  valley  between  Estevan  and  Roche  Percee*  Sir* 
James  Hector,  travelling  with  the  Palliser  Expedition  in 
1857,  measured  a  section  in  the  Souris  River  valley  and 


assigned  a 

Tertiary 

age  to 

t  he 

lignite 

coal  measures* 

Underground 

mining  of 

the  lignite 

b  egan 

in  1887 

after 

geological 

invest iga ti 

nns  in 

the 

area  by  A 

•R*C*  Selwyn  in 

1873  and  G.M 

•  Dawson  in 

1875.  A 

c  omprehensi ve 

description  of 

the  exposure 

s  along  the 

va lley 

walls 

of  the 

Souris 

Rive  r , 

Long  Creek 

and  Short 

Creek 

is 

found  in 

Dowling  ( 

1904  >. 

MacLean  (  1918  )  added  fu 

rther  descriptions  of 

the  coal 

seams 
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in  the  Souris  River  coalfield. 


Although  the  exploration  for  coal  was  most  active  in 
the  Estevan  area,  the  Tertiary  stratigraphy  of  southern 
Saskatchewan  was  deciphered  in  the  Cypress  Hills  area  where 
Upper  Cre taceous— Paleocene  transition  strata  are  well 

exposed.  The  term  Ravenscrag  was  first  applied  to  coal- 
bearing  strata  at  Ravenscrag  Butte  by  Davis  (1918).  Fraser 
<  1935)  recognized  two  parts  of  the  Ravenscrag  of 
differing  aspect  and  age.  The  term  Upper  Ravenscrag  was  used 
to  describe  the  Paleocene  coal— bearing  part  while  Lower 
Ravenscrag  was  applied  to  the  Cretaceous  non— coal— bearing 
part.  Fraser  ej^  (  1935  )  mapped  the  distribution  of  the 

Ravenscrag  Formation  in  southern  Saskatchewan  and  correlated 
the  Isolated  Souris  River  sections  as  part  of  their  Upper 


Ravenscrag  Formation.  Their  regional  geology 


and 


stratigraphy  until  recently  has  been  the  only  source  of 
information  on  the  formation.  In  1935,  Berry  described  the 
flora  of  the  Whitemud  and  Ravenscrag  Formations,  and 
Furnival  (  1946  )  renamed  the  Lower  Ravenscrag  the  Frenchman 
Formation  on  the  basis  of  his  work  in  the  Cypress  Hills 


area . 


More  recently,  Holter  (1972)  made  extensive  use  of  coal 
exploration  boreholes  to  correlate  individual  coal  zones  in 
the  Estevan  coalfield.  The  petrographic  composition  of  some 
seams  in  that  area  was  studied  by  Cameron  and  Birmingham 


(1971)  and  Broughton  (1972).  Whitaker  (1967,  1974a,  1974b) 


ugnj 

*- 
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produced  a  series  of  maps  and  sections  illustrating  the 
general  stratigraphy  and  structure  of  southern  Saskatchewan 
for  the  purpose  of  groundwater  exploration*  This  thesis  is 
based  on  data  collected  and  interpreted  during  the  lignite 
evaluation  program  described  by  Whitaker  (1972)* 

The  first  coal  resource  estimates  in  southern 
Saskatchewan  were  produced  by  Dowling  (  1913)*  Mackay  (  1947), 
using  more  realistic  criteria,  produced  new  estimates  for 
the  Royal  Commission  on  Coal*  Additional  knowledge  of  the 
area  and  revised  resource  classification  schemes  reduced 
Mackay1 s  estimates  by  one— half  to  those  of  Latour  and 
Chrismas  (1970)*  Finally,  on  the  basis  of  700  boreholes 
drilled  in  1972  and  1973,  and  500  boreholes  of  comparable 
quality,  Broughton  e  %  a l *  (  1974)  estimated  the  total 
quantity  of  measured  and  inferred  coal  resources  of 
immediate  economic  interest  to  be  5*7  x  109  tons* 

2*2  fiQgjoPftl  Qgglogy 

The  areas  which  form  the  basis  for  this  study  are 
located  in  the  Interior  Plains  physiographic  province  and 
within  the  Western  Canada  sedimentary  basin*  They  lie  on  the 
northern  flank  of  the  Williston  Basin,  an  ellipsoidal 
structural  depression  centered  in  North  Dakota*  Christopher 
et  a l *  (1973)  summarized  the  the  geological  setting  of  the 
area*  The  Precambrian  erosional  surface  dips  southwe stward 


at  12  to  20  feet  per  mile  as  an  irregular  monocline,  except 


' 


I 
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in  the  southeast  where  dips  to  the  south  of  up  to  50  feet 
per  mile  are  recorded.  Local  domal  features  originated  with 
or  were  rejuvenated  by  Laramlde  movements  (Kent  and  Simpson* 
1973  ). 

A  basal  clastic  sequence  records  marine  transgressions 
across  the  craton  in  Late  Cambrian  and  Early  Ordovician 
time.  Subsequent  to  the  development  of  the  Williston  Basin 
in  Early  Ordovician  time*  a  thick  carbonate— evaporite 
sequence  was  deposited  (Webb*  1964).  Although  interrupted  by 
an  hiatus  between  the  Silurian  and  Middle  Devonian* 
sedimentation  continued  Into  the  Mi ssi ssi pi an*  to  be 
followed  by  a  long  period  of  uplift  and  erosion. 


During  the  Jurassic*  evaporites*  carbonates*  shales  and 
sandstones  were  deposited  in  the  area.  Following  uplift  in 
the  Cordillera  in  Late  Jurassic/early  Cretaceous  time,  large 
volumes  of  elastics  were  deposited  around  the  margins  of  the 
Cretaceous  Boreal  and  Gulfian  seas.  Non— marine  sandstones 
and  marine  shales  intertongued  as  the  seas  transgressed, 
joined  and  finally  receded  to  the  south  and  east  as  further 
late  Cretaceous  and  Tertiary  uplift  and  vulcanism  occurred 
to  the  west. 


The  early  Paleocene  Cannonball  Sea 
represents  the  last  marine  incursion  into 
plains  before  the  area  was  covered  by  an  ext 


( Lemke 
the 

ensi ve 


,  1960) 
interior 
alluvial 


plain  represented  by  the  Ravenscrag  Formation.  As  a  result 
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of  the  late  Tertiary  L&ramide  orogeny*  southern  Saskatchewan 
was  uplifted  and  subjected  to  erosion*  Thin  stream— bed 
conglomerates  represented  by  the  Cypress  Hills  and  Wood 
Mountain  Formations  accumulated  with  erosional  unconformity 
on  the  Ravenscrag  Formation  (Taylor*  1964)*  Sands  and 
gravels  of  the  Empress  Group  (Whitaker  and  Christiansen* 
1972)  occur  as  valley  fill  deposits  in  pre— glacial  and  early 
glacial  valleys  on  the  bedrock  surface*  An  extensive  blanket 
of  glacial  drift*  locally  in  excess  of  900  feet  thick*  now 
covers  much  of  the  area* 

2*3  geplflgy  aJL  .th.s  R&xsknssjz&a  Eargali-qn 

Furnival  (  1946)  and  Fraser  e t  al *  (  1935)  provide  the 
most  recent  published  descriptions  of  the  Ravenscrag  in 
southern  Saskatchewan* 

The  Ravenscrag  Formation  consists  of  grey  and  buff- 
weathering  shales*  fine  sandstones*  siltstones*  and  sandy 
shales  and  conformably  overlies  the  late  Upper  Cretaceous 
Frenchman  Formation  (Furnival*  1946)*  Lignite  seams  are 
abundant  and  thin  laminae  of  carbonized  plant  material 
occur*  Fossil  tap  roots*  lithlfied  tree  trunks  and  stumps  in 
growth  position  have  also  been  found  in  the  formation* 
Bentonite  beds  form  a  minor  component  of  the  succession;  the 
sandstones  are  fine  to  very  fine  grained  and  exhibit  cross 


bedding  and  local  channeling  and  filling* 


. 

. 

. 


17 


Fossils  collected  Iron  the  Savenscrag  include  fresh* 
water  invertebra te,  fish*  turtle  and  crocodile  remains* 
Berry  (1935)  concluded  that  the  floral  assemblages  were 
representative  of  a  humid*  temperate  climate*  and  Fraser  e t 
&l.  (  1935)  Interpreted  the  deposl tional  environment  as  an 
alluvial  plain  covered  by  forests*  swamps  and  marshes* 
Azimuth  measurements  on  cross-stratification  (Carrigy,  1971) 
and  heavy-mineral  analyses  ( Kahmani  and  Lerbekmo*  1975)  have 
suggested  sediment  transport  from  the  northwest  for  the 
westernmost  exposures  of  the  Ravenscrag* 


2.4  gjBfll  Ggg.logy 


Thick  seams  of  lignite  (greater  than  3  feet)  in  the 
Ravenscrag  are  confined  to  four  major  basins  (Broughton  ^  i; 
al«  «  1974)*  Three  to  six  seams  occur  near  the  surface  in 
each  of  the  basins  and  are  separated  by  clastic  intervals 
varying  from  25  to  150  feet  thick*  All  samples  analysed  by 
Tibbetts  (1975)  are  lignite* 


UacLean  (1917*  p*  38)  characterized  seams  in  the  Souris 
River  valley  as? 

'••••not  continuous  through  the  area*  In  some 
places  they  split  up»  or  pinch  out  completely, 
to  start  again  at  some  distance  from  the  last 
observed  occurence*" 

Dowling  (  1904,  p*  9)  correlated  groups  of  seams  he  termed 
horizons  which: 


"•••may  contain  several  seams  but  in  each* 
workable  seams  occur  in  places." 


. 

. 


. 
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Throughout  an  area  of  coal  deposition)  individual  coal  seams 
may  thicken,  thin,  split  into  thinner  units  separated  by 
rock  or  disappear  altogether*  The  term  "zone"  is  used  here 
to  refer  to  any  interval  of  strata  which  contains  coal  in 
one  or  more  seams  and  which  can  be  correlated  over  a  wider 
area  than  any  individual  seam* 


Petrographic  studies  have  been  restricted  to  the 
Estevan  zone  in  the  Estevan  coalfield*  Cameron  and 
Birmingham  (1971)  recognized  ten  petrographic  units  varying 
from  3  to  12  inches  in  thickness  in  a  column  of  lignite  98 
inches  in  height*  They  concluded  that  the  maceral 
composition  favoured  a  dominantly  forest-moor  environment  of 
deposition  with  indications  of  reed— moor  and  open-water 
environments*  Broughton  (1972)  recognized  seven 
environmental  transitions  in  a  145  inch  column  of  the 
Estevan  zone*  Six  transitions  occured  from  the  dominant 
forest— moor  facies  to  the  less  frequent  reed— moor  and  one 
occured  from  the  reed— moor  to  the  open— water  facies* 


The  general  character  of  the  Ravenscrag  Formation  is 
very  similar  to  coal  measures  in  the  Sydney  Coalfield 
described  by  Haequebard  and  Donaldson  (1964)*  Their 
conceptual  model  of  a  flood— plain  depositional  environment 
is  an  appropriate  framework  for  interpreting  coal 
exploration  data  for  the  study  areas* 


In  Haequebard  and  Donaldson's 


model,  river  channels 


. 
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meandered  across  broad  areas  of  low  relief  carrying  sediment 
eroded  from  an  uplifted  area*  Peat  accumulated  in  large 
forested  swamps  and  marshes  which  covered  the  flood— plain* 
Overbank  fines  and  sandy  deposited  in  crevasse  splays  as  a 
result  of  floods  interrupted  peat  accumulation*  Occasionally 
channel  migration  or  switching  resulted  in  erosion  of  narrow 
channels  in  a  peat  deposity  and  during  pulses  of  orogenic 
uplifty  floods  of  elastics  burled  and  preserved  the  peat 
accumulations*  As  the  sediment  flow  decreasedy  plants 
stabilized  the  channels  and  peat  accumulation  resumed* 

2*5  Ealfl  Cgffpllfl.tl.9n 

All  exploration  data  used  in  this  study  are  from 
boreholes*  Lithologies  penetrated  were  described  by  a 
geologist  and  most  holes  were  logged  at  a  scale  of  l,,=20l 
using  single— point  resistance,  spontaneous  potentialy  gamma 
ray,  gamma— gamma  density  and  caliper  logs*  Gamma  ray  and 
gamma— gamma  density  logs  at  an  expanded  scale  of  1"  to  2' 
were  run  over  the  major  seam  intervals* 

A  coal  zone  comprises  deposi tional  events  which  can  be 
identified  and  correlated  using  subsurface  techniques* 
Letters  were  assigned  to  the  identifiable  beds  of  clean  coal 
in  a  zone  in  alphabetic  order  from  the  topmost  bed  down*  For 
operational  consistency,  each  bed  of  coal  one  foot  thick  or 
thicker  was  identified  as  a  seamy  with  seams  separated  by 
barren  strata  at  least  two  feet  thick*  Inorganic  partings 


\ 
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less  than  two  feet  in  thickness  often  occur  in  seams  and  the 
combination  of  coal  beds  in  the  seam  was  represented  by  the 
range  of  letters*  The  Estevan  C— G  seam  would  contain  the 
coal  beds  C»  D,  Ey  F  and  G*  Figure  3  illustrates  this 
nomenclature  as  applied  to  a  coal  zone* 

The  thickness  of  a  seam  as  recorded  is  the  total 
interval  from  the  top  to  the  base  of  the  seam*  The  net  coal 
thicknessy  in  contrasty  is  the  seam  thickness  minus  any  non- 
combustible  material  which  could  occur  in  beds  less  than  2 
feet  thick*  These  terms  are  illustrated  in  Figure  4*  Data  on 
the  location  and  identity  of  boreholes  and  the  stratigraphy 
and  coal  seams  penetrated  were  codedy  keypunched  and  entered 
into  a  computer  data  file  which  formed  the  basis  for  this 
study* 

Because  of  thickness  and  depth  co nsi derat i ons y  not  all 
seams  in  a  zone  at  a  borehole  might  reasonably  be  mined* 
Therefore  a  method  was  was  devised  to  combine  selected  seams 
of  a  zone  into  a  single  unit  termed  a  potential  mining  zone* 
Seams  which  were  mathematically  acceptable  at  each  borehole 
according  to  the  criteria  described  below  were  treated  as  a 
unit  having  a  bottom  depth  equal  to  the  base  of  the  lowest 
coal  seam  chosen  and  a  net  coal  thickness  equal  to  the  total 
of  the  net  coal  thicknesses  of  the  accepted  seams*  Criteria 
used  for  selection  of  seams  from  each  zone  for  inclusion  in 
the  potential  mining  zone  in  each  borehole  were: 

1*  Seams  with  five  feet  or  more  of  net  coal  are 
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FIGURE  3.  SEAM  NOMENCLATURE  APPLIED  TO  A  HYPOTHETICAL  COAL  ZONE 


22 


ZONE 

THICKNESS 
(in  feet) 

28 

28 

SEAM 

THICKNESS 
(in  feet) 

n 

Ci 

Ci 

o 

NET 

COAL 
(in  feet) 

Ci 

ci 

in 

SECTION 
(in  feet) 

Ci 

00 

mm 

r 

m 

::: 

if 

Ci 

00 

ci 

1 

.  ‘  • 

TOTALS 

LITHOLOGY 

:  :  COO! 

Coal 

I  v  pm 

o  a 
o  o 

U  U 

o 

_C 

to 

"a 

o 

U 

NAME 

MODIFIER 

< 

:  CO  tj:  ; 

Q  UJ 

O 

ZONE 

NAME 

NVA31S3 

SEAM 

NUMBER 

- 

V:.'::  ";:  'CN  . 

Figure  4.  Diagrammatic  explanation  of  coal  thickness  terms. 
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always  Included* 

2*  Seams  with  less  than  five  feet  of  net  coal  but 
with  three  feet  or  more  are  included  if  the 
incremental  overburden  ratio  (overburden  or 
rock  material  above  the  seam  and  below  the 
next  higher  seam  or  surface  divided  by  the 
coal  thickness)  is  less  than  152 1,  or  if 
underlain  by  a  seam  fitting  criterion  1* 

3*  Seams  with  less  than  three  feet  of  net  coal 

and  with  two  feet  or  more  are  included  when 

underlain  by  a  seam  fitting  either  criteria  1 
or  2 • 

OS 

When  overlain  by  an  acceptable  seam  provided 
that  the  rock  interval  between  the  seams  is 
less  than  twice  the  thickness  of  the  lower 
seam  • 

The  validity  of  the  above  criteria  is  dependent  upon 

assumptions  about  the  value  and  end  use  of  the  coal  and  the 
present ly—prac ticable  mining  methods  employed  in  southern 
Saskatchewan • 


In  some  borehole 
under  the  above  sele 
erosion  or  non— depos 
boreholes t  the  entire  a 
missing  was  excluded 


s  no 

seams  in  a 
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c  t  ion 

criteria* 

Where 

evidence  of 

1 1  ion 

c  ou Id  be 

found 

in  adjacent 

rea  in 

which  the 

zone  was 

too  thin  or 

from 

consideration  in 

creating  an 

' 


24 


isopach  map  and  in  the  subsequent  resource  quantity 
calculation*  Where  the  lack  of  acceptable  seams  appeared  to 
be  restricted  to  a  single  borehole  and  did  not  extend  to  the 
surrounding  areay  a  zero  value  for  net  coal  in  the  potential 
mining  zone  was  used  for  resource  quantity  estimation  and 
for  isopaching  the  net  coal  values* 

Most  coal  zones  are  truncated  by  the  bedrock  erosional 
surface  which  developed  after  the  Paleocene  and  prior  to  the 
onset  of  glaciation*  To  predict  the  location  of  the  subcrop 
of  the  potential  mining  zone  beneath  glacial  deposits^ 
structure  contour  maps  of  the  top  of  the  zone  and  the 
bedrock  surface  were  produced  from  borehole  data  by  computer 
programs  using  various  algorithms  for  estimating  the  value 
of  a  surface  at  the  nodes  of  a  regular  grid*  Although  grid 
intersections  form  a  regular  geographic  basis  for  surface 
comparisons^  the  choice  of  a  grid  spacing  involved  a  trade¬ 
off  between  the  desired  level  of  detail  in  a  contour  map  and 
the  expense  of  creating  very  fine  grids*  A  grid  spacing  of 
500  metres  was  selected  for  the  Estevan  area  whereas  a  100 
metre  distance  between  nodes  was  used  for  the  small-scale 
area*  The  grids  were  then  contoured  and  algebraically 
compared  using  the  SURFACE  II  graphics  package* 

Few  exploration  boreholes  were  drilled  in  stream 
valleys  where  recent  erosion  has  removed  the  overlying  till 
and  exposed  the  bedrock  surface  and  often  the  coal  zones 
proper*  Bedrock  surface  grids  produced  from  borehole 
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information 

thus 

usually 

did 

not 

reflect 

such 

valleys* 

Selected  contour 

lines 

from 

the 

publi shed 

1:50 

000  scale 

topographic 

maps 

were  digitized 

and 

a  surface 

topography 

grid  was  produced*  This  surface  topography  grid  was  compared 
to  the  bedrock  surface  grid,  and  where  a  surface  elevation 

value  was  lower  than  the  corresponding  estimated  bedrock 

/ 

elevation,  the  surface  elevation  was  substituted  into  the 
bedrock  grid,  adding  the  detailed  information  from  the 
surface  topography  to  the  bedrock  grid  data*  The  location  of 
the  subcrop  of  the  coal  zone  on  the  bedrock  surface  was  then 
estimated  by  algebraically  comparing  the  modified  bedrock 
grid  and  the  coal  zone  structure  grid* 

2*6  £gtexaP  A-Cgft 

The  Estevan  coalfield  is  located  in  southeastern 
Saskatchewan  (Figure  2);  the  location  and  distribution  of 
the  340  boreholes  in  the  coalfield  used  in  this  study  are 
shown  in  Figure  5* 

The  region  is  flat  to  gently  rolling  prairie  dissected 
by  the  Souris  River,  Long  Creek  and  Short  Creek*  The  naxlmur 
relief  from  the  bottom  of  the  river  valleys  to  the 
International  Boundary  in  the  southeast  is  about  200  feet 
(67  m)  (Figure  6)*  Meneley  .aJL*  (  1957)  believe  that  the 
Souris  River  and  Long  Creek  occupy  meltwater  channels  which 
were  marginal  to  the  Weyburn  glacial  lobe  at  two  different 
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The  Estevan  coalfield  underlies  a  plateau  of  bedrock 
bounded  to  the  west,  north  and  east  by  the  valley  cut  by  the 
northeast— flowing  preglacial  Missouri  River  ( Whitaker r 
1974a;  Meneley  e  t  a\*  *  J957)*  The  elevation  of  the  bedrock 
plateau  varies  from  1825  feet  to  1950  feet  above  sea  level 
(Figure  7)*  The  northwest— trending  channel  of  a  preglacial 
river  tributary  to  the  major  river  isolates  a  bedrock  high 
to  the  southwest  of  the  main  upland*  Sands  and  gravels  of 


the  Empress  Group  cover  the  bottoms  of  most  of 


the 


preglacial  channels*  Incision  of  the  seltwater  channel  now 
occupied  by  the  Souris  River  exposed  bedrock  (including 
coal)  to  a  depth  of  100  feet  or  more  below  the  upland  level 
in  an  east— west  band  about  a  mile  in  width* 


At  least  five  mineable  coal  zones  have  been  identified 
in  the  Estevan  coalfield  (Hotter,  1972)*  Because  of  their 
wide  areal  extent,  contrasting  seam  character  and 
availability  of  data,  the  Estevan  zone  and  the  underlying 
Boundary  zone  were  chosen  for  study* 


2*7  Eale yap  z ane 


The  Estevan  zone  was  studied  in  the  area  bounded  on  the 
south  by  the  International  Boundary,  on  the  east  by  the 
67C  000  metre  UTM  easting  and  to  the  north  and  west  by 
estimated  subcrop  boundaries*  The  zone  was  intersected  in 
176  boreholes  (shown  in  Figure  8  )J  incomplete  sections  of 
the  Estevan  zone  are  preserved  in  some  15  boreholes  located 
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near  or 

the  zone 
In  these 
used  for 


outside  the  subcrop  edge  where  the  upper  seams  in 
have  been  removed  by  glacial  or  preglacial  erosion* 
latter  boreholes  the  top  of  the  remaining  zone  was 
structural  interpretation* 


A  sand  body,  interpreted  as  a  stream  channel 
contemporaneous  with  coal  accumulation*  was  mapped  for  about 
9  km  along  the  north  boundary  of  Township  1*  Range  7  W2 
(Figure  8)*  After  most  of  the  zone  was  deposited,  the 
channel  was  apparently  abandoned  and  filled  with  two  upper 
coal  beds*  The  area  in  which  the  coal  filling  was  too  thin 
to  be  mineable  or  is  missing  was  considered  as  an  area  of 
non-deposition  and  was  excluded  from  resource  quantity 
calculations* 


The  zone  surface  generally  strikes  east— nor theast  and 
dips  southeasterly  at  20  feet  per  mile  with  the  dip 
steepening  to  greater  than  30  feet  per  mile  to  the  southwest 
along  the  east  flank  of  a  sout h^plunging  trough  underlying 
Roche  Percee  and  Bienfait  (Figure  9)*  The  dip  of  the  zone 
appears  to  flatten  out  to  12  feet  per  mile  or  less  at  the 
subcrop  under  and  southwest  of  Estevan*  A  small  south- 
plunging  anticline  is  present  in  the  southwest  corner  of  the 
map  southwest  of  the  bedrock  channel*  The  anticline  may 
continue  northeast  across  the  bedrock  channel  although  seam 
correlations  are  not  sufficiently  reliable  to  confirm  this* 


Local  structure  is  present  as  northwest— trending  swales 
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of  low  amplitude  having  an  apparent  wavelength  of  one  to  two 


m 


lies.  The  orientation  of  these  structural  features 


perpendicular  to  the  direction  of  Pleistocene  ice  movement 
and  near  seam  subcrop  and  valley  walls  suggests  a  glacial 
origin • 


The  Estevan  zone  is  composed  of  up  to  9  coal  beds  which 
can  be  correlated  separately  or  combined  in  seams*  An 
average  of  two  seams  is  found  in  each  borehole  although  as 
many  as  four  may  occur*  The  average  thickness  of  the 
complete  zone  is  14  feet*  and  the  net  coal  thickness  varies 
from  2  feet  (the  minimum  acceptable  according  to  the 


c  r 


iteria)  to  18*4  feet  with  an  average  of  9*5  feet*  An 


average  of  1*4  seams  per  borehole  were  accepted  in 
calculating  the  total  net  coal  thicknesses*  The  net  coal 
thickness  in  the  zone  tends  to  be  greater  in  the  west  and 


no  r 


th  (Figure  10),  but  each  area  contains  exceptions  to  this 


;eneralization* 


2.8  Boungarz  Zan?. 

The  Boundary  coal  zone  lies  20  to  86  (average  60)  feet 
below  the  Estevan  zone*  Bounded  on  the  south  by  the 
International  Boundary,  and  to  the  north  and  west  by 
subcrop,  the  Boundary  zone  was  limited  by  sand  depostion  to 
the  east*  These  limits  and  the  distribution  of  the  160 
boreholes  which  intersect  the  zone  are  shown  on  Figure  11* 
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Due  to  "the  wide  spacing  of  boreholes  over  "the  southeast 
part  of  the  areay  structure  contours  on  the  top  of  the 
Boundary  zone  (Figure  12)  lack  the  detail  of  the  local 
structure  found  in  the  Estevan  zone*  The  zone  generally 
strikes  east— northeast  and  dips  to  the  south  at  25  feet  per 
mile  but  the  dip  becomes  more  easterly  on  the  west  side  of  a 
slight  depression  located  between  Bienfait  and  Roche  Percee* 
The  local  features  on  the  regional  structure  echo  the 
northwest- trending  swales  on  the  overlying  Estevan  zone  but 
appear  to  be  more  subdued  and  are  smaller  in  dimension* 

Up  to  twelve  coal  beds  in  the  zone  have  been 
tentatively  correlated  among  boreholes*  The  zone  comprises 
one  to  five  seams  with  an  average  of  two  seams  in  each 
boreholey  and  a  total  zone  thickness  of  about  9  feet* 
Thickness  of  the  potential  mining  zone  was  calculated  in  151 
hotes  and  set  to  zero  in  9  holes  in  which  the  seams  were  too 
thin*  Up  to  3  seams  were  selected  in  some  boreholes  but  the 
average  number  of  seams  ,used  in  forming  the  potential  mining 
zone  was  1*2  per  borehole* 

Zero  values  of  net  coal  thickness  in  several  boreholes 
are  apparently  due  to  local  thinning  of  seams  below  the 
cutoff  value  for  the  selection  criteria  and  could  not 
clearly  be  interpreted  as  a  change  in  lithology  as  in  the 
case  of  the  Estevan  zone*  These  zero  net  coal  values  were 
Included  in  the  construction  of  the  isopach  map  and 
subsequent  resource  estimation*  Net  coal  thickness  in  the 
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potential  mining  zone  varies  from  zero  to  a  maximum  of  13 
feet  with  an  average  of  5#3  feet# 

Thicker  sections  of  coal  occur  in  the  west  and  north 
but  the  thickness  distribution  as  illustrated  by  the  isopach 
of  net  coal  (Figure  13)  appears  to  be  very  irregular#  The 
numerous  small  closures  of  the  isopach  lines  is  in  part  due 
to  the  Interpolation  algorithm  used  in  the  SURFACE  II 
contouring  package# 


2.9  Ssaie  A-Ega 


Data  in  the  form  of  over  420  borehole  logs  were 
available  for  the  study  of  small  scale  variability  of  a  coal 
zone  in  an  area  about  three  miles  square*  Figure  14 
illustrates  the  distribution  of  the  control  boreholes  and 
surveyed  section  lines# 

The  area  lies  on  the  gently  sloping  south  wall  of  a 
broad  valley#  The  surface  topography  is  shown  in  Figure  15# 
A  misfit  stream  which  flows  along  the  thalweg  of  fhe  valley 
across  the  northeast  corner  of  the  area  is  Joined  by  two 
intermittent  streams  from  the  south#  The  bedrock  surface 
(Figure  16)  is  composed  of  two  elements*  a  broad*  flat 
valley  floor  to  the  northeast  bounded  to  the  southwest  by  a 
steeper  valley  wall#  The  break  in  slope  of  the  bedrock 
surface  trends  southeast  across  the  centre  of  the  area#  The 
valley  wall  was  dissected  by  two  tributary  streams  draining 
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INCOMPLETE  SECTION  OF  COAL  ZONE 

COAL  ZONE  IS  tllSSINO 


FIGURE  14. 


DISTRIBUTION  OF  BOREHOLES 
IN  THE  SHALL  SCALE  AREA. 
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FIGURE  15.  SURFACE  TOPOGRAPHY  OF  THE 

SMALL  SCALE  AREA. 
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FIGURE  16.  BEDROCK  TOPOGRAPHY  OF  THE 

SMALL  SCALE  AREA. 
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the  upland  to  the  south*  Glacial  till  and  discontinuous 
sands  and  gravels  of  the  Empress  Group  overlie  the  bedrock 
surface  and  range  in  thickness  from  70  feet  in  the  valley  to 
20  feet  on  the  valley  wall* 

The  coal  zone  has  a  regional  dip  of  30  feet  per  mile  to 
the  northeast  but  the  structure  is  locally  complex  (Figure 
17  )•  A  small  depression  40  to  60  feet  deep  occurs  in  the 
southeast  corner  of  the  area  and  appears  to  be  bounded  to 
the  west  and  north  by  a  north-trending  anticlinal  flexure* 
Local  normal  faults  with  offsets  of  up  to  20  feet  were 
inferred  along  the  trend  of  the  structure  but  could  not  be 
traced  over  distances  of  one- eighth  of  a  mile*  A  small 
northeast-plunging  synclinal  structure  crosses  the  centre  of 
the  area* 

Similar  structures  have  been  ascribed  to  salt  collapse 
(DeMllle,  1964;  Christiansen,  1967)*  Solution  and  removal  of 
salt  in  the  Prairie  Evaporite  Formation  of  Devonian  age  has 
resulted  in  subsidence  of  the  overlying  sediments  and 
structural  adjustment  along  normal  faults*  Salt  solution  of 
the  Prairie  Evaporite  was  active  during  the  Tertiary  and  may 
have  controlled  sediment  dispersal  and  coal  accumulation 
locally  through  subsidence  of  negative  elements* 

The  coal  zone  is  composed  of  up  to  six  beds  but  the 
lower  two  are  thin  and  erratic  in  distribution*  A  complete 


section  of  the  zone  consists  of  an  upper  seam  separated  by  a 


' 


■ 


. 

■ 


44 


FIGURE  17.  STRUCTURE  CONTOURS  ON  THE 

BASE  OF  THE  COAL  ZONE. 
SMALL  SCALE  AREA. 
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sand  parting  from 
thickness  locally  in 
southeastf  a  fact  which  may 
positive  element  at  the  time  of 
salt  which  dissolved  after  the 
The  distribution  of  the  zone 
in  Figure  14* 


parting  increases  in 
depression  in  the 
that  the  area  was  a 
deposition*  a  remnant  of 
deposition  of  the  Ravenscrag. 
in  the  boreholes  is  indicated 


a  lower  seam*  The 
the  structural 
indicate 


Figure  18  is  an  interpretation  of  the  areal 
distribution  of  the  coal  zone*  A  sand-filled  channel  which 
existed  during  coal  deposition  has  been  has  been  inferred 
from  boreholes  to  occur  in  a  narrow  east— west  bahd*  The  zone 
has  been  truncated  by  preglacial  erosion  at  the  foot  of  the 
valley  wall  except  where  it  was  preserved  in  the  previously- 
mentioned  synclinal  flexure.  Downcutting  in  tributary 
valleys  removed  the  coat  from  two  isolated  areas* 


Net  coal  thickness  in  the  zon 
affected  by  erosion  and  non— deposi tion 
from  less  than  8  feet  in  the  southwest 
in  the  central  area,  with  an  average  of 
thickenings  occur  adjacent  to  the  sand 
appear  to  be  erratic  over  the  remainder 


ey  excluding 
( Figure  19  )  * 
to  more  than 
about  12  fee 
channel  but 
of  the  area* 


va lues 
ranges 
16  feet 
t*  Local 
va lues 
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FIGURE  18. 


SAND  CHANNEL  AND  EROSIONAL 
AREAS  IN  THE  COAL  ZONE. 
SMALL  SCALE  AREA. 
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FIGURE  19.  I SGPflCH  OF  NET  COAL  THICKNESS 

OF  THE  COAL  ZONE  IN  THE 
SMALL  SCALE  AREA. 
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CHAPTER  3 


3*  Slali 

In  1965f  Koch  end  Gomez  employed  statistical  techniques 
to  estimate  the  thickness  of  some  lignite  beds  in  Texas  but 
they  did  not  consider  the  problem  of  area  estimation*  As 
results i  they  presented  the  arithmetic  mean  of  the  thickness 
observations  in  boreholes  and  confidence  intervals  for  the 
mean  based  on  the  t  distribution  and  the  sample  variance* 

For  their  results  to  be  meaningful)  four  assumptions 
must  be  satisfied  (Koch  and  Gomez*  1965*  p  4)2 

( 1  )  each  seam  is  homogeneous 

(2)  the  boreholes  are  drilled  at  random 

(3)  results  from  different  boreholes  are  statistically 
1 ndependent 

(4)  the  observations  are  normally  distributed* 

The  third  assumption  is  the  most  significant  and  most 
likely  to  be  violated  in  reality*  For  example*  consider 
measurements  of  the  height  of  adult  males  passing  a  street 
corner*  The  height  of  the  next  man  to  pass  the  corner  cannot 
be  predicted  by  the  value  of  a  particular  measurement  taken 
previously*  and  is  statistically  independent*  There  is  no  £ 
priori  reason  to  suspect  that  the  next  individual  will  be 
tall  or  short  but  rather  that  the  height  of  all  adult  males 
forms  a  frequency  distribution  characterized  by  certain 
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parameters  such  as  a  mean  and  variance#  According  to  the 
third  assumption  seam  thickness  as  observed  in  boreholes  is 
similar  in  character  to  the  above  process.  To  continue  the 
analogy(  consider  a  field  in  which  men  have  been  assembled 
and  are  now  standing  shoulder  to  shoulder.  An  isoline  map  of 
the  height  of  the  men  would  appear  chaotic.  GiantSf  midgets 
and  those  of  average  height  would  be  randomly  arranged  in 
the  field. 

This  model  does  not  fit  previous  observations  of  coal 
seam  thicknesses  or  geological  models  of  peat  deposition.  In 
outcrops  along  valleys  and  in  minest  coal  seam  thicknesses 
are  observed  to  vary  relatively  continuously.  If  considered 
as  a  mineable  zone  in  which  the  contribution  of  several 
seams  may  be  combined)  the  net  coal  thickness  has  a  regional 
continuity  with  local  thickenings  and  thinnings  of  coal.  A 
simplistic  analogy  of  a  similar  process  is  men  passing  by 
the  street  corner  in  groups  of  similar  heights  the  relative 
position  at  which  an  individual  arrives  in  the  sequence  of 
groups  is  significant  for  accurate  prediction  of  the 
individual's  height  as  location  is  significant  for  accurate 
prediction  of  coal  seam  thickness. 

The  processes  of  peat  accumulation  at  any  one  point  are 
controlled  by  many  factors.  A  coal  zone  probably  represents 
ten  or  more  discrete  configurations  of  peatland  and  alluvial 
channels  and  the  coal  thickness  at  any  point  represents  the 
sum  of  these  processes.  As  each  of  the  major  peat 


. 


50 


deposi tional  events  occurred  (and  shifted)  reasonably 
uniformly  over  an  area,  the  sum  of  these  processes  has  a 
relatively  continuous  spatial  structure.  Minor  variations  in 
bed  thickness  due  to  differential  compaction,  growth 
patterns  of  vegetation,  erosional  channelling  and  water 
table  level  are  impossible  to  detect  in  boreholes  spaced  on 
an  average  of  more  than  half  a  mile  apart*  Any  model  used  to 
describe  coal  thickness  variation  must  take  into  account  not 
only  the  seeming  random  component  but  also  a  structured 
aspect • 

From  the  above  discussion  it  follows  that  not  all  the 
results  from  different  boreholes  are  statistically 
independent*  Observations  taken  close  to  one  another  are 
more  likely  to  be  similar  than  observations  separated  by  a 
great  distance.  Statistical  dependence  might  be  expected 
between  thickness  measurements  from  closely  spaced  boreholes 
while  observations  widely  spaced  would  be  statistically 
independent*  The  dependence  or  correlation  between  samples 
of  the  same  phenomenon  at  different  locations  is  called 
spatial  autocorrelation  (Agterberg,  1974,  p*  314)* 


The  assumption  that 
(Assumption  2)  implies 
as  seam  thickness  that 
observations  is  a  biased 
Areas  which  are  more 


the  boreholes  are  drilled  at  random 


for  an  autoco rrelated  variable  such 


the  arithmetic  mean  of  the 


estimator  of  the  average  thickness, 
heavily  sampled  than  the  average 

i  s 


density  contribute  more  weight  to  the  mean 


than 


. 
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appropriate  for  an  unbiased  estimator* 

To  sunmarizei  the  assumptions  made  by  Koch  and  Gomez 
(1965)  are  not  intuitively  appealing*  The  coal  seam 
thickness  is  not  randomly  distributed  spatially*  and  there 
are  priori  reasons  to  expect  patterns  to  occur*  On  the 
other  hand*  the  coal  thickness  cannot  be  predicted  exactly 
because  variaticns  occur  at  random  locations  due  to  small 
scale  processes*  Finally*  an  irregular  sampling  pattern 
results  in  a  fragmented  picture  of  the  structured  and  random 
aspects  of  the  phenomenon* 

Equations  will  be  introduced  in  the  remainder  of  the 
chapter  to  describe  the  properties  of  some  random  processes 
and  the  estimation  of  their  parameters*  In  these  formulae* 
•Z'  will  refer  to  a  structural  surface  elevation,  net  coal 
thickness  or  any  other  continuously  mappable  variable  which 
exhibits  spatial  correlation.  The  location  of  any  point  in 
the  plane  will  be  referenced  as  ux"  and  that  of  the  i— th 
point  location  as  "x( i  )"•  The  i-th  measurement  has  a  value 
"Z( i  ) " •  The  location  of  the  point  or  area  which  is  to  be 
estimated  is  Mx( o )  and  "Z(o)M  is  the  value  of  the  estimate* 


3 •  1  Stationary  Stochastic 


Mathematical  models  of  processes  which  could  be  used  to 
describe  coal  thickness  variation  over  an  area  are 
stationary  stochastic  processes  in  the  plane  (Matern,  1960, 
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p.  11).  A  stochastic  process  Is  stationary  if  the  expected 
mean  value  is  a  constant  and  the  covariance  function  of  two 
observations  depends  on  the  distance  between  the  points  and 
not  their  location* 

This  may  be  described  by  the  following  equations.  Let 
Z( x )  be  a  random  variable  at  the  location  x  in  the  plane. 
Then  the  mean  function  of  Z  is 

m(  x  )  =  E[Z(x)]  =  constant  (3.1) 

and  the  covariance  function  is  defined  as 

C(x,y)  =  E[  Z(  x  )Z(  y  )  ]-m(  x  )m(  y  )  (3.2) 

where  E  is  used  to  denote  •expectation*.  The  symbols  x  and  y 
represent  point  locations. 

The  autocorrelation  function  of  the  process  is 

r(  x  *  y  )  =  C(x,  y  )•[  C(  x,  x  )C(  y,  y  )  ]“  l/2 

which  can  be  rewritten  as 

r(  x  t  y  )  =  C(x.y)  <3.3  ) 

cr 

where  O’2  is  the  population  variance. 

Since  the  covariance  function  depends  only  on  the  distance  h 

between  x  and  y»  (3.3)  can  be  rewritten  as 

r(h)  =  C( h)  =  QLhl  where  h  =  x-y 

c(  o  )  <r 

C< o  )  is  the  covariance  of  an  observation  to  itself  and  it  is 
theoretically  the  population  standard  deviat ion, CT • 

This  model  is  in  wide  use  in  time  series  analysis  for 
measurements  taken  at  regular  time  intervals.  The  distance 


. 
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between  measurements  is  called  the  lag*  The  behaviour  of  the 
autocorrelation  function  with  distance  can  be  studied  using 
a  correlogram •  An  experimental  autocorrelation  coefficient 
r( k  )  is  calculated  for  all  pairs  of  measurements  k  lags 
apart  and  plotted  on  the  ordinate  with  the  lag  as  the 
abscissa*  In  general*  the  coefficients  decrease  from  1  to  0 
and  even  attain  negative  values  as  lag  increases*  Examples 
of  corre log rams  of  geologic  variables  may  be  found  in 
Agterberg  (1970;  1974,  Chapter  10)  and  Krige  &l.  (1969). 

Permissible  models  of  autocorrelation  functions  in  two 
dimensions  have  been  discussed  by  Matern  (1960),  Krige  et 
Al*  (1969)  and  Agterberg  (1970)*  Agterberg  concluded  that 
the  form 

r(h)  =  Cexp(-ah)  (3*4) 
could  be  fitted  to  experimental  correlograms  as  a  model  of 
the  autocorrelation  function*  If  the  coefficient  C  is  less 
than  one*  the  phenomenon  is  discontinuous  near  the  origin* 
The  coefficient  a  is  a  constant  which  scales  the  decrease  in 
correlation  to  the  distance  h* 

Matern  (1960)  reviewed  previous  work  on  spatial 
variation  and  studied  the  efficiency  of  unrestricted  random 
sampling,  stratified  random  sampling,  Latin  square  designs 
and  systematic  sampling  of  stationary  and  isotropic 
processes  with  a  decreasing  correlation  function*  He  applied 
his  results  to  optimum  sampling  of  forest  tracts* 
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Ag terberg  and  Chung  (  1973)  computed  a  two-dimensional 
anist ropic  autocorrelation  function  for  "the  sulfur  content 
in  the  Harbour  seam*  Sydney  coalfield*  They  used  the 
experimental  autocorrelation  coeff icierts  to  predict  the 
sulfur  content  using  the  best  linear  unbiased  estimator 
(BLUE)*  The  BLUE  of  a  statistical  parameter  is  a  linear 
combination  of  all  the  known  values,  Z(i);  it  is  unbiased; 
it  has  the  minimum  variance  of  all  possible  linear,  unbiased 
estimators • 

In  the  case  of  a  stationary  random  function*  the  BLUE 
of  an  unknown  point  Z(  o  )  may  be  characterized  by 

Z(o>  =  >^a(i)Z(i)  i  =  1»N  (3.5) 

where  ^  ^a(  J  )C(  ij  )—  p  =  C(oi) 

^a(  i  )  =  1  (3.6) 

the  a( i  )  are  the  weights  of  the  Z( i  )  used  to  form  the 
weighted  average  Z(o).  C(ij)  is  the  covariance  between 
sample  i  and  sample  j.  C(oi)  is  the  covariance  between  the 
point  to  be  predicted  and  sample  i«  The  Lagrangian 
multiplier*  p  is  used  to  obtain  a  minimum.  The  second 
condition,  that  the  sum  of  the  weights  is  unity,  assures 
unbiasedness . 

If  the  mean  of  a  variable  is  non— stationary  or  only 
weakly  statlonaty*  the  correlogram  has  several  operational 
drawbacks  as  a  measure  of  spatial  variation.  The  variance  of 
a  non— s ta tl onary  process  increases  with  the  size  of  the 


field  sampled  and  if  the  variance 


is  not  finite* 


the 


. 
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covariance  and  the  autocorrelation  function  are  undefined* 

3.2  Intrinsic  SftndSEl  Functions 

A  weaker  hypothesis  to  replace  the  assumption  of 
stationarity  is  the  intrinsic  hypothesis  (Matheron,  1971,  p. 
53).  He  defines  an  intrinsic  random  function  as  having 
increments  which  are  wide  sense  stationary*  That  is# 

E[  Z(  x+h  )-Z(  x  )  ]  =m(h)  (3*7) 

and 

E[(  Z(  x+h  )-Z(x  )  )2  ]  =  27(h)  (3.8) 
where  E  denotes  the  expectation.  The  function  7(h)  on  the 
right  hand  side  of  3.8  is  called  the  semi vari og ram  or  more 
commonly  the  variogram.  The  expression  m(h)  is  the  expected 
difference  between  observations  a  distance  h  apart  and  is 
called  the  drift. 

This  formulation  appears  to  be  more  convenient  In  that 
it  can  accomodate  a  linear  trend  in  Z(x)  and  7(h)  is  defined 
even  If  no  a  priori  variance  exists  for  the  variable.  If  the 
variograa  remains  bounded  with  distance  (finite  variance) 
the  relationship 

7(h)  =  C(  o  )— C(  h  )  (3.9) 
exists  between  the  variogram  and  covariance.  The  phenomenon 
is  wide-sense  stationary  if  the  drift  is  constant  and  the 
variance  is  finite.  The  variogram  (serial  variation  function 
of  Jowett  (1955))  is  a  basic  tool  in  Matheron* s  Theory  of 
Regionalized  Variables. 


It  is  used  to  characterize  the 
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spatial  variation 
as  very  irregular 
aspect  and  a 
measured  in  one  or 
variables • 


of  regionalized  variables  which  he  defines 
functions  of  location  composed  of  a  random 
structured  aspect*  Continuous  phenomena 
more  spatial  dimensions  are  regionalized 


Using  the  relationship  between  the  variogram  and 
covariance  (3*9)*  the  BLUE  of  a  stationary  random  function 


C  3. 10  ) 
( 3.11  ) 


(3*6)  may  be  rewritten  as 

Z(  o  )  =  y\(  i  )Z(  i  ) 

^ji(  j  )T(  i  j  )+p  =  T(  oi  )  i  =  ltN 

2>‘ 1 » =  1 

where  *y(  i  j  )  is  the  value  of  the  variogram  between  the  pair 
of  samples  i  and  J  and  T( oi  )  is  similarly  the  value  of  the 
variogram  for  the  point  to  be  estimated  o  and  the  sample 
point  i*  This  method  was  termed  kriglng  by  Matheron  in 
honour  of  D.G*  Krige  who  used  a  similar  form  of  regression 
of  the  actual  grade  of  ore  blocks  on  the  grade  predicted 
from  samples*  Although  Matheron* s  theory  is  verj  rich  in 
applicat ions »  the  following  discussion  will  be  confined  to 
two  topics!  structural  analysis  and  kriglng* 


3*3  structural  Analysis 


The  structural  properties 
using  experimental  variograms* 
(3*8)  can  never  be  obtained 
taken  over  a  finite  area  may  be 


of  a  variable  may 
The  theoretical 
but  the  variogram 
computed  using  th 


be  studied 
variogram 
of  samples 


formula 
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y(  h  )  =  1/2N  [  Z(  x(  1  )  )-Z(  x(  i+h  )  )  ]2  (3.12) 
N  is  "the  number  of  palps  of  samples  separated  by  the  vector 
h  with  a  direction  8  where 

tan0  =  y/x  (3.13) 

and 

h2  _  x2+y2  (3.14) 

If  sample  points  are  regularly  distributed  along  a  line 
of  length  L  then  the  experimental  variogram  can  be  estimated 
for  h  up  to  L/2  but  is  most  reliable  near  the  origin.  When 
sample  points  are  irregularly  spaced  in  the  plane*  pairs  of 
points  are  classed  into  distance  intervals*  h±dh.  The  class 
width*  2dh  is  selected  so  that  most  classes  have  30  pairs  of 
observations  or  more  to  ensure  a  stable  variogram  estimate. 
For  each  distance  class*  the  drift 

m(h)  =  1/n7>(  x(  1  )  )-Z(  x(  i  )+h  )  ]*  (3.15) 
the  variogram*  and  the  average  distance  between  pairs  are 
calculated.  If  the  variogram  is  estimated  without  regard  for 
the  direction  of  the  vector  h*  it  is  called  the  average  or 
isotropic  variogram. 

Anisotropy  in  the  spatial  structure  of  a  variable  can 
be  detected  if  variograms  are  computed  in  several 
directions.  Since  sample  points  are  rarely  available  on 
lines  in  the  direction  of  interest,  a  pair  is  included  if 
the  vector  is  within  an  angular  tolerance  of  the  desired 
direction.  An  angular  tolerance  of  ±30°  is  acceptable  since 
only  general  trends  are  of  interest  in  the  structural 
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analysis*  Figure  20  illustrates  how  pairs  of  points  are 
included  in-  experimental  variograms* 

The  purpose  of  structural  analysis  is  to  fit  a 
geologically  reasonable  model  of  the  variograra  function  to 
the  observed  experimental  va riograms( s  )  •  The  model  fitted  is 
used  for  the  purpose  of  estimation  and  although  many 
functions  might  fit  the  experimental  points*  Matheron  (1971) 
suggests  a  restricted  range  derived  from  theoretical 
processes*  Two  classes  of  models  are  used  in  practice! 
variograms  which  are  bounded  and  variograms  which  continue 
to  increase  »ith  distance  (Figure  21)* 

Bounded  variograms  rise  to  a  sill  and  remain  constant 
with  increasing  distance*  They  are  characteristic  of 
phenomena  which  have  a  finite  variance  (the  sill)»  are  wide- 
sense  stationary  and  obey  the  relation  (3*9)*  The  distance 
at  which  the  model  rises  to  the  sill  is  the  range  or  zone  of 
influence  of  the  process.  Points  separated  by  a  distance 
greater  than  the  range  are  statistically  independent*  If  the 
range  measured  in  one  direction  differs  markedly  from  that 
observed  in  another  direction,  anistropy  in  the  variable  is 
indicated* 

One  of  the  most  common  bounded  models  is  the  spherical 
scheme  defined  as! 

r(h)*cr^lx-l  llll  h  <  a  (3*16) 

2  a  2  a3J 


r(  h ) 


c 


h  >  a 
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DISTANCE  AND  ANGULAR  SELECTION  CRITERIA  FOR  PAIRS  OF  POINTS 
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MODELS  ANISOTROPY 
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NESTED  STRUCTURES 


Figure  21. Characteristics  of  Variograms 
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where  a  is  the  range  of  the  variogram  and  C  is  the  variance 
or  sill* 

Serra  (1968)  has  shown  that  such  transitional  processes 
can  be  observed  at  the  scale  of  microns,  centimetres,  metres 
and  kilometres  in  the  grade  of  iron  in  the  Lorraine  Iron 
Formation*  The  nature  of  processes  with  ranges  shorter  than 
the  closest  sampling  interval  cannot  be  measured  on  an 
experimental  variogram  but  can  be  inferred  in  the  model 
fitted  to  the  data*  An  experimental  variogram  which  clearly 
shows  a  break  in  slope  may  be  modelled  as  two  or  more  nested 
or  intermeshed  structures  of  different  ranges  and 

intensities  (Figure  21)* 

Variograms  which  are  not  bounded  characterize  non¬ 
stationary  phenomena  with  a  non— constant  drift  (3*15)  and  a 
variance  that  is  a  function  of  the  size  of  the  area  sampled. 
Such  variograms  are  often  linear  and  may  indicate  ani stropy 
by  a  different  slope  in  each  direction*  Because  the 

variogram  is  not  bounded,  there  is  no  range  or  zone  of 
influence  of  sample  points*  Transition  phenomena  with  ranges 
shorter  than  the  sampling  interval  appear  in  experimental 
variograms  as  a  non— zero  intercept  even  though  .such 
processes  are  spatially  continuous*  A  linear  variogram  is 
modelled  using  the  equation 

y(h)  =  bh  (3.17) 

where  b  is  the  slope  and  h  is  the  distance* 


' 
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The  behaviour  of  the  variogram  near  the  origin  is  at 
the  same  time  important  and  difficult  to  determine*  The 
interpretation  has  a  critical-  effect  on  the  weighting  given 
to  samples  in  linear  prediction  and  on  the  value  of  the 
estimation  variance*  The  detailed  information  required  is 
available  only  where  continuous  exposures  of  the  geological 
phenomena  are  available*  The  dense  drilling  required  to 
collect  such  information  for  buried  coal  deposits  is  often 
not  justified*  Modelling  the  nature  of  the  variogram  between 
the  origin  and  the  first  few  experimental  points  is 
necessarily  an  art  of  some  skill  and  subjectivity* 

Figure  21  indicates  the  shape  of  four  types  of 
variograms  near  the  origin*  The  continuous  type  is 
representative  of  bed  thickness  and  surface  elevations 
measured  without  error*  Errors  in  measurement  or  slight 
discontinuities  in  the  phenomena  introduce  what  is  referred 
to  as  a  nugget  effect*  The  random  type  of  variogram 
corresponds  to  the  spatial  structure  of  the  height  of  adult 
males  randomly  arranged  in  a  field  or  the  thickness  of  a 
coal  seam  under  the  assumptions  made  by  Koch  and  Gomez 
(1965)  and  discussed  earlier* 


The  nugget  effect  is  so  called  because  of  the  effect 
that  very  concentrated  nuggets*  small  relative  to  the  volume 
of  a  sample*  have  on  the  experimental  variogram  of  gold 
values*  Analogous  problems  occur  in  coal  quality 
observations  which  are  dependent  upon  the  physical  volume  of 
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the  sample •  Sabourln  (1975)  has  noted  a  nugget  effect  in  the 
sulfur  content  of  coal  zones  in  Cape  Breton* 


Since  borehole  observations  are  point  observations  of 
continuous  surfacesi  the  experimental  variograms  were 
interpreted  in  this  study  as  continuous  near  the  origin  with 
the  possibility  of  a  nugget  effect  from  errors  in 
observation*  If  a  nugget  effect  is  present,  the  models 
previously  described  (3*16  and  3*17)  become 

Y(h)  =  Co+CitiJ  h  -  1  il33  h  <  a  (3.18) 

2  a  2  a3 

Y(h)  =  C0 +Ct  h  >  a 


and 


Y(h)  =  CQ+bh 


(3.19) 


where  Co 
variance 


is  the  value  of  the  nugget  effect  and  Co+Ci  is 
for  the  spherical  model  (3*18). 


the 


The  purpose  of  structural  analysis  in  this  application 
is  to  identify  variogram  models  which  can  be  used  for 
solving  estimation  problems.  Although  structural  analysis  is 
a  powerful  tool  for  data  analysis,  caution  must  be  exercised 
in  the  geological  interpretation  of  experimental  variograms. 
Dutta  and  Rao  (1977)  discussed  the  problems  encountered  in 
the  structural  analysis  of  mixed  populations.  David  (1974, 
pp.  68-142)  provides  an  excellent  discussion  of  practical 
techniques  of  variography  as  well  as  the  inherent  pitfalls. 
Some  points  on  the  experimental  variogram  may  be  Ignored, 


and  it  can  be  used  for  estimation,  if  the  model  fitted  is 
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consistent  with  most  of  the  data* 

3*4  srigiaa 

Once  a  model  of  the  variogram  function  has  been  chosen 
by  structural  analysis*  the  model  can  be  used  to  calculate 
the  estimation  variance  for  various  estimators*  The 
estimation  variance  is  the  variance  of  the  distribution  of 
possible  errors  made  in  estimating  a  true  unknown  value  Z  by 
a  sample  value  or  a  combination  of  sample  values*  Z** 
Hopefully*  an  estimator  will  be  unbiased  such  that 
E[ z-z*  ]  =  0 

and  the  estimation  variance 
E[ ( Z— Z*  )2  ] 
will  be  a  minimum* 

The  estimation  variance  of  extending  the  value  of  a 
point  sample  to  a  line  and  to  regular  geometrical  shapes  has 
been  charted  for  standard  variogram  models  (Matheron*  1971, 
pp.  87-95).  FORTRAN  subroutines  used  to  evaluate  these 
auxiliary  functions  for  the  spherical  scheme  (equation  3*16) 
have  been  published  by  Clark  (19761*  Some  commonly  used 
estimators  are: 

(1)  the  value  of  a  central  sample  extended  to  the  enclosing 
rectangle • 

(  2  )  the  mean  of  the  corner  samples  extended  to  the 
rectangle • 

(3)  the  arithmetic  mean  of  samples  taken  on  a  rectangular 


■ 


' 

- 


X  o 

* 


65 


gr  id  • 

(4)  the  arithmetic  mean  of  samples  collected  by  random 
stratified  sampling  with  rectangular  strata. 


The  estimation  variances  of  more  complex  geometrical 
relationships  between  the  samples  and  the  area  estimated  can 
be  calculated  as  combinations  of  the  elementary  shapes. 
Kriging  has  the  property  of  being  the  estimator  which  has 
the  minimum  estimation  variance  of  all  the  possible  linear 
combinations  of  sample  points  and  is  therefore  preferable  to 
most  of  the  other  estimators  available. 


Local  estimation  is  the  prediction  of  the  unknown  value 
of  Z  at  a  point  x( o )  or  the  average  value  Z  over  a  regular 
polygon  which  is  small  in  comparison  to  the  area  sampled. 
The  kriging  system  in  this  case  is: 

Z(o)  =  £a(  i  )Z(i  ) 

£]a(  j  )Y(  i  J  1+  M  =  T(oi)  i=l,N 

i>  =  i 


The  quantity  T(  i J  )  is  the  value  of  the  variogram  between 
sample  points  i  and  J  and  is  easily  evaluated  from  the 
variogram  model.  The  a(  1  )  are  the  weights  and  /z  is  the 
Lagrangian  multiplier.  When  x( o  )  is  a  point,  Y(oi)  is  the 
value  of  the  variogram  between  the  sample  point  x( i  )  and  the 
point  to  be  predicted  x(o).  When  x( o )  is  a  regular  polygon, 
y( oi  )  becomes  Y(oi),  the  mean  value  of  the  variogram  between 
points  in  the  polygon  and  the  sample  point  x(i).  *y(  oi  )  can 
be  evaluated  using  the  auxiliary  functions  discussed 
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earlier. 


The  system  of  equations  is  composed  of  a  matrix  of 


sample  to  sample  variogram  values 


S  = 


Til  >12  >13  •••  >ln  1 

>2  i  >22  >23  ...  T^n  1 


*>n  f  >n  2  Vn  3  ...  >nn  1 
111  10 


a  vector  of  sample  to  unknown  variogram  values 


U= 


>01 

>02 


Ton 


and  the  vector  of  weights  which  are  the  unknowns 


A= 


al 
a  2 

• 

an 


The  system  is  written  as 
SA=U 

and  solved  for  A  by  rewriting  as 
A=S~ 1 U 

where  S— 1  is  the  inverse  of  the  matrix  S.  The  vector  A  is 
then  used  to  weight  the  samples  vector  Z  to  form  an  estimate 
of  Z(o).  The  estimation  variance  by  kriging  is 

(T2-Xla(  1  )Y(  oi  °  ^ 

where  V( o  )  is  the  mean  variogram  value  of  all  points  in  the 
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polygon  to  be  estimated  to  all  other  points  in  the  polygon# 
If  x(  o  )  is  a  point, 

Y(  o  )=0 
ot herwi se 


Auxiliary  functions  can  be  used  to  evaluate  Y( o )  for  regular 
shapes. 

Global  estimation  is  the  prediction  of  the  average 
value  Z  of  a  variable  over  the  areal  extent  of  that 
variable.  The  kriging  system  of  equations  is  identical  to 
that  in  the  local  estimation  problem.  The  values  of  Y  are 
estimated  by  the  sum 


where  x( j  )  is  a  point  on  a  square  grid  located  at  random 
over  the  area  and  l  is  the  number  of  grid  points  within  the 
polygon  x(o).  Since  the  boundaries  of  most  areas  are 
irregular ,  the  grid  spacing  must  be  such  that  there  are 
greater  than  200  grid  intersections  in  an  area  for  a  stable 
estimate  of  Y( oi  ).  Y( o  )  is  estimated  by  a  similar  discrete 
summation  method. 

3.5  Nog— Stationary  Fh£flgJ8£J3a 

In  contrast  to  a  wide-sense  stationary  proeess  which 
has  a  constant  mean  function  and  the  weakly  stationary 
process  which  has  a  constant  mean  for  a  local  area,  the  mean 
function  of  some  variables  clearly  depends  upon  the 
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location.  Structural  surfaces  of  inclined  and  disturbed 
strata  are  examples  of  variables  where 

m  ( x  )  =  E[f(x)]  is  not  constant. 

Agterberg  (1970),  Parsley  (1971)  and  Matheron  (1971)  have 
suggested  models  of  such  non~stationary  phenomena. 

Agterberg  (1970)  proposed  that  each  observation  of  a 
geologic  attribute  in  space  is  the  sum  of  three  components: 
DATA=TREND+SIGNAL+NOISE 

Variation  called  "signal"  at  a  large  scale  might  become 
"trend"  at  a  small  scale  and  similarly  "noise"  might  become 
"signal"  as  more  detailed  sampling  became  available.  The 
"trend"  component  is  the  large  scale  "regional"  variations 
which  may  be  represented  by  the  fitting  of  deterministic 
polynomials  to  the  data.  The  signal,  according  to  Agterberg 
(1970,  p.  120)  is  "...weakly  stat ionary. • •  [and] 

.. .characterized  by  its  autocorrelation  function".  Noise  is 
the  purely  random  and  uncorrelated  component  caused  by  rapid 
variation  over  short  distances  and  measurement  errors.  His 
approach  to  modelling  spatial  variation,  as  evidenced  by  the 
terminology,  is  adapted  from  time-series  analysis.  The  use 
of  deterministic  polynomials  to  represent  the  trend  is  the 
familiar  trend  surface  analysis  (Whitten,  1975).  The 
statistical  model  for  the  signal  is  a  weakly  stationary 
random  function  of  distance. 

Parsley  (1971)  described  the  components  of  variation  in 
mapped  geologic  data  in  the  following  manner: 
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TOTAL=REGIONAL+LOCAL+RES IDUAL 

To  separate  the  three  components,  he  employed  a  low  order 
trend  surface  to  estimate  the  regional  component  and  then 
fitted  successively  higher  order  polynomials  to  deviations 
from  the  regional  trend  until  serial-correlation 
coefficients  computed  for  residuals  indicated  that  the 
residuals  were  not  au tocor relat ed •  The  use  of  a  high— order 
trend  surface  to  model  the  local  variation  contrasts  with 
Agterberg's  choice  of  a  stochastic  or  random  function# 


Matheron's  geostatist ical  model  of  non-s tat  ionary 
processes  has  been  described  by  Delf iner  ( 1976,  p#  51 )  as 
DATA  =  DRIFT  +  RESIDUAL 

The  drift  is  his  term  for  the  mean  function  and  the 
fluctuation  or  residual  is  the  random  component#  The  drift 
is  similar  to  the  trend  component  of  Agterberg's  model  and 
is  represented  by  local  deterministic  polynomials#  The 
residual  is  assumed  to  be  a  random  function  with  a  spatial 
structure  which  can  be  modelled  by  a  variogram  function# 


In  practice,  the  variogram  can  be  a  biased  estimator  of 
the  covariance  function  of  the  residual  component  because 
the  drift  m(  x  )  is  not  known  but  must  be  estimated  (Delfiner, 
1976).  If  the  drift  is  evident  as  in  the  homoclinical  dip  of 
a  bed,  the  experimental  variogram  of  the  residuals  along  the 
strike  is  probably  a  good  estimator#  If  trends  are  present 
but  are  not  constant  across  the  area  studied,  a  more 


comprehensive  method  for  identifying 


spatial  structure  is 
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based  on  the  theory  of  intrinsic  random  functions  of  order  k 
(IRF—k)  and  their  generalized  covariances  (Mathercn,  1973)* 

Delfiner  (1976)  provides  an  introduction  to  linear 
prediction  using  IRF— k  theory*  An  intrinsic  random  function 
of  order  k  is  a  random  function  whose  kth  order  increments 
are  weakly  stationary*  The  spatial  structure  of  the  kth 
order  increments  are  estimated  as  isotropic  generalized 
covariances  appropriate  lor  the  order  k*  Once  an  order  k  and 
the  coefficients  of  the  generalized  covariance  have  been 
chosen  ,  a  modified  form  of  krlging  called  universal  kriging 
is  used  for  local  estimation*  Universality  conditions  are 
equations  of  polynomials  up  to  order  k  representing  the 
drift  effect  which  are  added  to  the  basic  kriging  system* 
The  coefficients  of  the  polynomial  terms  and  the  estimation 
variance  can  be  computed  as  part  of  the  solution* 

Structural  analysis  of  non-s tati onary  phenomena  using 
the  variogram  tends  to  be  a  biased  estimator  except  in  rare 
situations*  An  automated  procedure  has  been  designed  for 
selecting  the  order  k  and  the  parameters  of  the  generalized 
covariance*  The  method  consists  of  estimating  each  sample 
point  by  all  other  sample  points  assuming  orders  0,  1  or  2 
and  comparing  the  mean  squared  error  of  each  order*  The 
order  with  the  lowest  mean  squared  error  is  the  best 
predictor*  Once  the  order  k  is  selected,  the  generalized 
covariance  model  which  most  closely  approximates  the 
observed  estimation  variance  is  chosen*  The  resulting 
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es  time  to  r  is  empirically  "the  best  estimator  for  the 

observations • 

If  samples  are  taken  from  heterogeneous  surfaces 
composed  of  rough  and  smooth  areas*  the  resulting  estimated 
generalized  covariance  is  suitable  for  neither  area  and  may 
be  biased  if  one  of  the  areas  is  more  heavily  sampled* 
Furthermore*  the  model  selected  is  only  one  of  many  models 
which  fit  the  data*  Another  model  may  be  more  appealing  when 
other  information  and  conceptual  geological  models  are  taken 
into  account* 

3«7  Areal  Methgdg 

The  preceding  sections  have  been  concerned  with  methods 
of  estimating  the  value  of  a  continuous  variable  at  a 
location  or  the  average  of  a  continuous  variable  over  an 
area*  This  section  examines  methods  of  estimating  the  area 
underlain  by  coal* 

One  method  is  to  estimate  the  areal  proportion  of  coal 
underlying  the  map  area  by 

p  =  n/N 

where  n  is  the  number  of  boreholes  intersecting  coal  and  N 
is  the  total  number  of  boreholes*  If  p  is  distributed  as  the 
binomial  distribution,  confidence  intervals  for  the  estimate 
can  be  made*  This  model  is  based  on  the  naive  assumption 
that  observations  are  statistically  Independent  —  ie*  that 
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an  observation  located  a  short  distance  (say  5  feet)  away 
from  a  previous  observation  of  coal  has  the  probability  p  of 
being  coal.  Switzer  (1976)  suggested  that  the  estimator 
could  be  improved  by  systematic  sampling  on  a  square  grid 
and  that  the  estimation  variance  could  be  calculated  from  a 
form  of  autocorrelation  function  adapted  from  Matern  (I960). 
His  method  requires  a  systematic  sample  pattern  which  is 
rarely  available  in  practice. 


Matheron  (1971)  has  proposed  another  estimator  of  area 
based  on  systematic  sampling.  A  square  cell  is  centred 
around  each  sample  point  and  the  area  estimate  is 

S=Na 

where  N  is  the  number  of  cells  in  which  the  borehole 
intersected  coal  and  a  is  the  area  of  a  cell.  Since  the 
location  of  the  boundary  is  not  known*  the  estimator  has  a 
relative  variance  dependent  on  the  length  of  the  boundary 
relative  to  the  area.  This  is  approximated  for  the  isotropic 
case  by  the  number  of  boundary  segments  between  cells 
containing  coal  and  cells  without  coel.  For  the  method  to  be 
useful f  the  boundary  of  the  area  must  be  outlined  by 
sampling  and  the  sample  spacing  must  be  reasonable  for  the 
expected  variability  of  the  phenomenon.  The  formula  for 
estimating  the  relative  variance  of  the  area  estimator  is 

Nt  >  N2  ( 3. 20  ) 


=  —I 

S2  N2 


1  N2+0.061ULU2 

6  N2 


This  method  is  once  again  dependent  upon  the  availability  of 


systematic  sampling. 
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Neither  of  these  methods  use  the  available  structural 
data*  In  the  present  study t  much  of  the  area  underlain  by 
coal  is  bounded  at  subcrop  so  that  the  prediction  of 
location  of  the  subcrop  is  a  form  of  area  estimation*  Grid 
manipulation  methods  permit  surface  comparison  to  facilitate 
prediction  of  the  subcrop  line*  The  subcrop  line  is  the  zero 
contour  on  an  isopach  of  the  depth  of  a  seam  below  the 
bedrock  surface*  The  isopach  can  be  represented  as  the 
difference 

fiedrock  —  .Seam 

for  each  grid  node*  The  variance  of  the  possible  error  made 
in  estimating  a  surface  at  a  grid  node  can  be  produced  by 
universal  kriging*  The  variance  of  the  sum  of  two 
independent  random  variables  is  the  sum  of  the  individual 
variances 

o-2  CB+S)  =  <r  2  (  B  )  +<r2(S) 

Since  the  bedrock  surface  and  the  seam  structure  were  formed 
by  different  processes  and  are  estimated  from  different 
control  points  using  different  models  they  are  considered 
here  as  independent*  Confidence  intervals  about  the  isopach 
value  at  each  grid  intersection  can  be  determined  if  the 
error  of  estimation  is  assumed  to  be  normally  distributed 
with  a  zero  mean  and  a  variance  equal  to  the  estimation 
variance.  The  location  of  the  subcrop  line  for  each 
confidence  level  can  be  displayed  on  a  contour  map*  The  area 
can  then  be  planimetered  or  integrated  digitally  and  a 
distribution  of  probability  versus  area  can  be  computed* 


. 
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Such  a  distribution  would  not  be  regular  in  shape  since  the 
area  estimates  may  skew  quite  strongly* 

The  method  of  area  estimation  using  subcrop  lines  is 
only  a  partial  solution  to  the  total  problem  since  in  some 
cases  the  areal  extent  of  coal  is  bounded  by  a  lithologic 
change  to  clastic  facies*  All  three  coal  zones  described 
previously  have  lithologic  boundaries  interpreted  where  they 
grade  into  sand  bodies  or  other  lithologies*  The 
interpretation  of  the  nature  and  extent  of  such  boundaries 
is  a  subjective  judgement  made  by  a  geologist  based  on  often 
limited  borehole  data*  Statistical  inference  about  the 
location  of  such  boundaries  based  on  objective  probability 
is  impossible*  The  objective  probability  of  an  event  may 
only  be  determined  either  theoretically  from  the  physical 
characteristics  of  a  mathematically-defined  process  or 
empirically  as  the  result  of  a  large  number  of  trials* 

Subjective  probability  offers  a  solution  to  statistical 
inference  when  the  objective  probability  approach  cannot  be 
used*  Subjective  probability  is  the  degree  of  belief  or 
personal  opinion  of  the  likelihood  that  an  event  will  occur. 
Expert  geological  opinion  can  be  used  to  define  a  cumulative 
probability  distribution  of  the  form  P(  X  >  x)  for  many  types 
of  parameters*  For  the  areal  extent  of  a  coal  seam*  there 
are  only  three  parameters  of  the  distribution  which  can  be 
reasonably  defined:  the  maximum*  the  minimum  and  the  mode  or 
most  likely  value*  The  triangular  distribution  (Newendorp, 
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1975»  p.  273)  can  be  used  to  create  a  continuous  cumulative 
probability  distribution  based  on  these  three  parameters* 


•  <. 

. 
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LL  Estevaa  Area 


The  Location  of  coal  exploration  boreholes  commonly 
does  not  form  a  systematic  pattern»  a  fact  which  precludes 
the  use  of  statistical  techniques  based  on  a  regular 
geometrical  configuration  of  samples.  Non-sys temat ic 
sampling  of  spatially  correlated  variables  gives  rise  to  the 
problem  of  weighting  observations  to  form  an  unbiased 
estimate.  Kriging  and  subjective  probability  appear  to  offer 
solutions  to  the  problem  of  providing  quantitative 
uncertainty  limits  to  coal  quantity  estimates  from 
irregularly  distributed  borehole  data. 


4 .  1  Location  gflal  Sflfc><?£flJ2 

In  section  3.6*  a  method  was  proposed  which  would 
predict  the  location  of  coal  subcrop  lines  and  allow 
confidence  limits  to  be  placed  on  the  estimated  locations.  A 
procedure  for  implementing  the  method  i s  as  follows: 

Cl)  Perform  a  structural  analysis  of  the  borehole  data  for 
the  bedrock  and  coal  zone  elevations.  Select  models  of 
the  experimental  variograms  and  the  order  k  if  a  drift 
is  indicated. 

(2)  Estimate  the  elevations  of  the  bedrock  and  the  coal  zone 
using  kriging  at  the  nodes  of  a  regular  grid.  Save  the 
estimation  variance  for  each  estimate  on  separate  grids. 
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(3)  Modify  the  bedrock  grid  to  reflect  stream  erosion  in 
recent  valleys.  This  is  performed  with  a  surface 
topography  grid  based  on  contours  digitized  from  maps. 

(4)  Subtract  the  coal  zone  surface  grid  from  the  modified 
bedrock  grid  to  form  a  grid  of  isopach  values.  Where 
grid  values  are  negative*  the  coal  zone  is  predicted  to 
lie  at  elevations  above  the  inferred  bedrock  surface  and 
is  assumed  to  have  been  removed  by  erosion.  The  zero 
contour  interpolated  from  the  isopach  grid  values  on  a 
contour  plot  is  interpreted  as  the  best  estimate  of  the 
subcrop  of  the  top  of  the  zone. 

(5)  Combine  the  estimation  variance  for  the  bedrock  surface 

with  the  estimation  variance  for  the  coal  zone  surface 
a t  each  grid  node  to  form  an  estimation  variance  for  the 
isopach  grid  value.  The  distribution  of  possible  error 
made  in  estimating  the  isopach  should  be  symmetrical 
about  a  zero  mean  value  (because  the  kriging  estimator 
is  unbiased)  with  a  variance  equal  to  the  estimation 
variance  (because  the  errors  of  estimation  for  the 
bedrock  and  the  coal  zone  are  assumed  to  be 

independent  )•  Although  there  have  been  few  studies  of 
actual  error  distributions  using  kriging,  David  (1974, 
p.  17)  suggests  that  the  normal  distribution  is  a  safe 
if  not  conservative  assumption. 

(6)  Subtract  (or  add)  one  or  more  standard  deviations  of  the 
estimation  error  at  each  grid  node  to  form  grids  of 
isopach  values,  each  grid  with  a  constant  cumulative 
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probability  that  the  real  isopach  thickness  at  a 
particular  node  is  greater  than  the  calculated  isopach 
value.  The  cumulative  probability  that  the  true  unknown 
value  (X)  is  greater  than  or  equal  to  a  value  (x)  is 
derived  from  a  cumulative  normal  distribution  with  a 
mean  equal  to  the  estimated  isopach  and  a  variance  equal 
to  the  estimation  variance  and  is  written 


p(x  2  x> 

For  a  distribution  with  a  mean  m  and  a  variance  s2  *  some 


values  of  P 

P(  X  2 

are: 

m-*2s  ) 

=  0.277 

P(X  2 

m— s  ) 

=  0.841 

P(X  2 

m  ) 

=  0.500 

P(X  2 

m+s  ) 

=  0.159 

Four  grids 

!  o  f 

isopach 

the  above  equations  and  are  contoured. 

The  probability  that  the  real  unknown  area  of  coal 
(X)  is  greater  than  the  area  enclosed  by  the  zero 
contour  line  (x)  is  the  same  as  the  probability  that  the 
real  isopach  thickness  at  a  grid  location  is  greater 
than  the  estimated  grid  value. 

(7)  Create  a  cumulative  probability  distribution  for  the 
area  enclosed  by  the  subcrop  lines.  Using  four  grids 
with  probability  values  P(  X  >  x)  =  .98*  .84»  .50  and 
•16*  the  areas  corresponding  to  the  subcrop  lines  can  be 


plotted  against  probability  to  form  a  distribution. 
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Structural  analysis  of  the  surfaces  was  based  on 
experimental  variograms  produced  by  a  program  VARIOGRAM  and 
an  automated  variogram- fit ting  routine  in  KRIGPACC,  a 
proprietary  software  system  of  Aquitaine  Company  of  Canada 
Ltd.  Five  experimental  variograms  were  calculated  for  each 
variable:  an  isotropic  variogram  and  four  directional 

variograms  made  at  45°  intervals.  The  analysis  performed  by 
KRIGPACC  reports  the  residual  sum  of  squares  for  kriging 
using  universality  conditions  of  order  0,  1  and  2.  This  can 

be  useful  to  select  an  order  of  polynomial  if  a  drift  is 
indicated.  For  each  order*  several  sets  of  coefficients  for 
the  generalized  covariance  are  suggested,  each  accompanied 
by  a  selection  criterion. 

Several  points  must  be  made  about  the  choice  of  an 
order  k  of  drift  and  a  model  of  the  variogram  of  residuals. 
First,  the  order  k  which  minimizes  the  residual  sum  of 
squares  is  not  always  the  wisest  choice.  The  minimum  number 
of  data  points  required  tc  estimate  a  point  using  kriging 
with  a  quadratic  drift  is  12  while  only  6  are  needed  if  a 
linear  drift  is  assumed  (Delfiner,  1975,  p.  65).  If  there  is 
only  a  marginal  difference  in  the  residual  sum  of  squares 
criteria  between  the  first  and  second  order  increments,  then 
selecting  the  linear  drift  would  permit  extrapolation  beyond 
the  limits  of  data  control. 

If  a  drift  is  indicated  in  the  analysis  and  selected, 
then  the  sample  weights  (and  hence  the  estimate)  do  not 
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depend  only  on  the  form  of  the  variogram  but  are  determined 
primarily  by  the  universality  conditions  and  hence  the 
geometry  of  the  data  points  and  the  point  to  be  estimated. 
The  choice  of  a  model  for  the  variogram  of  residuals  is  not 

critical  to  making  estimates  since 

"•••the  fundamental  inde term inabi li ty  will  only  affect 
the  variances  of  our  optimal  estimators  and  not  the 
estimators  themselves"  ( Matheron*  1971»  p*  189)* 

Care  was  taken  to  model  the  residual  variogram  as 

realistically  as  possible  since  the  variance  of  the 

estimator  was  to  be  used  in  constructing  confidence  limits. 

Experimental  variograms  for  the  elevation  of  the 

bedrock  surface  are  displayed  in  Figure  22.  The  average 
variogram  indicates  a  transition  type  phenomenon  rising  from 
the  origin  to  a  peak  of  approximately  2400  ft2  at  about  6  km 
and  dropping  down  to  1600  ft2  at  about  8  km.  The  cause  of 

the  "bump"  in  the  average  variogram  is  more  easily  seen  in 

the  directional  variograms.  The  variogram  oriented  at  45° 
(and  those  at  0°  and  90°  to  a  lesser  extent)  shows  the 
effect  of  a  systematic  drift  between  2  km  and  8  km.  The 
drift  is  caused  by  the  pairing  of  points  in  the  bedrock 
valley  southwest  of  Estevan  with  points  on  the  general  level 
of  the  bedrock  plateau.  The  variogram  calculated  at  135° 
parallel  to  the  preglacial  valley  shows  no  drift  effects  and 
has  beep  taken  as  the  variogram  of  the  residuals.  The 
residual  sum  of  squares  calculated  by  KRIGPACC  (Table  2) 
suggests  that  a  linear  drift  minimizes  the  estimation 
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DISTANCE  ,  METRES 

FIGURE  22.  BEDROCK  SURFACE,  VARIOGRAMS  AND  MODEL. 
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vari ance • 

The  kriging  parameters  lor  the  bedrock  surface  include: 
a  linear  drift  and  a  spherical  residual  variogram  with  no 
nugget  effect,  a  sill  of  1650  ft2  and  a  range  of  6500 
metres*  The  variogram  model  is  plotted  on  Figure  22*  Figure 
7  is  a  contour  map  of  the  bedrock  grid  estimated  by  the 
above  model  and  modified  for  the  surface  topography.  The 
standard  deviation  of  estimation  error  for  the  bedrock 
elevations  has  been  contoured  and  is  displayed  in  Figure  23* 
The  standard  deviation  ranges  from  a  theoretical  zero  at  the 
data  points  to  greater  than  30  feet  outside  the  drilled 
area*  The  distribution  of  values  is  controlled  by  the 
geometry  of  the  sample  points  and  the  data  void  located  in 
the  valley  of  the  Souris  River  is  marked  by  a  "ridge"  of 
high  standard  deviation. 

Experimental  variograms  for  the  top  of  the  Estevan  zone 
(Figure  24)  show  a  transition  phenomenon  with  no  nugget 
effect,  a  shoit  range  structure  with  a  sill  of  about  400  ft2 
at  800  metres  and  a  marked  anisotropy  for  distances  greater 
than  2500  metres*  The  directional  variogram  computed  along 
the  approximate  strike  of  the  surface  (  65a)  appears  hounded 
with  distance  whereas  variograms  in  ether  directions  show  a 
linear  increase  and  drift  components  beyond  a  distance  of 
2*5  km.  The  residual  sum  of  squares  (Table  2)  for  the  sample 
data  indicate  that  there  is  little  difference  between 
orders.  Because  a  systematic  trend  is  evident  in  the  dip  of 
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FIGURE  24.  ESTEVAN  ZONE  SURFACE,  VAR  I OGRAMS  AND  MODEL. 
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the  bed  and  insufficient  data  points  are  available  within  a 
2*5  km  neighbourhood  of  many  grid  nodes  in  the  southeast 
part  of  the  oapi  a  linear  drift  was  chosen*  The  variogram 
model  fitted  to  the  residuals  has  two  nested  spherical 
mode Is! 

Y(  h  )  =  T  i(  h  )+Y2(  h  ) 

and  the  parameters  of  each  (equation  3*18)  are 
Ci  =  225  ft2  ai  =  750  metres 

C2  =  1200  ft2  a2  =  15000  metres 

The  resulting  model.  is  plotted  with  the  experimental 
variograms  in  Figure  24*  Figure  9  is  the  contour  map  of  the 
Estevan  zone  structure  grid  as  estimated  by  universal 
kriging*  The  standard  deviation  of  the  estimation  error  for 
the  Estevan  zone  structure  is  shown  in  Figure  25*  The  areas 
of  poor  control  within  the  coal  field  have  a  standard 
deviation  of  20  to  25  feet*  In  the  area  drilled  at  centres 
one  mile  apart  north  of  the  river  and  east  of  Estevan*  the 
standard  deviation  varies  from  15  feet  to  20  feet  except  at 
locations  very  close  to  the  data  points*  The  standard  error 
of  the  Estevan  zone  surface  in  the  triangular  nose  of 
bedrock  south-west  of  the  bedrock  channel  varies  from  10 
feet  to  20  feet*  an  area  drilled  at  one— half  mile  centres* 

The  isopach  of  the  interval  between  the  bedrock  surface 
and  the  Estevan  zone  is  shown  in  Figure  26*  The  zero  contour 
line  is  the  best  estimate  of  the  location  of  the  subcrop  of 
the  coal  zone*  and  the  gradient  of  the  isopach  surface 
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FIGURE  25.  STANDARD  DEVIATION  OF  THE 

EST IN AT  ION  ERROR, 

ESTE  VAN  ZONE  SURFACE 
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indicates  the  angle  at  which  the  zone  is  truncated  at  the 
bedrock  unconformity#  This  angle  varies  from  3°  along  the 
margins  of  the  bedrock  channel  to  1°  at  the  erosional  edge 
to  the  northeast*  In  contrast)  the  estinated  surface  of  the 
Estevan  zone  is  almost  parallel  to  the  bedrock  unconformity 
under  the  city  of  Estevan  and  in  the  southeast  quarter  of 
Township  2,  Range  9*  The  standard  deviation  of  the 
estimation  error  for  the  isopach  is  displayed  in  Figure  27* 
The  values  range  from  zero  to  30  feet  with  the  most  common 
values  between  20  feet  to  25  feet* 

These  standard  deviations  were  used  to  create  lines 
representing  estimated  subcrop  locations  at  various 
probability  levels  (Figure  28)*  The  estimates  diverge 
markedly  in  the  area  to  the  northwest  but  are  closely  spaced 
and  sub— parallel  elsewhere*  Within  the  "best*'  subcrop  line 
to  the  northwest,  over  50  boreholes  drilled  at  close 
spacings  in  Township  2,  Ranges  8  and  9  did  not  intersect  the 
Estevan  zone*  Detailed  investigation  of  the  area  indicated 
that  the  bedrock  surface  is  horizontal  with  a  relief  of  10 
to  15  feet*  The  surface  of  the  Estevan  zone  which  can  be 
reconstructed  suffers  a  reversal  of  dip  from  the  regional 
trend  and  is  apparently  almost  horizontal*  The  structural 
interpretation  is  hampered  by  uncertainty  in  the  data 
arising  from  two  sources* 

In  shallow  boreholes,  the  section  studied  may  not 
include  any  recognizable  marker  horizons  and  the 
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correlations  made  are  very  subjective.  Additional  noise  may 
be  introduced  into  the  data  when  the  elevation  mapped  is  not 
a  true  stratigraphic  top  but  a  synthetic  horizon  calculated 
from  the  elevation  of  the  base  of  the  lowest  seam  accepted 
and  the  gross  thickness  of  coal  accepted  in  the  potential 
mining  zone. 

Clearly  the  "best"  subcrop  line  overestimated  the 
extent  of  the  Estevan  zone,  the  cause  of  this  being  the 
apparent  reversal  of  dip  in  the  southeast  corner  of  Township 
2»  Range  9  which  was  carried  over  to  the  data  void 
underlying  Estevan.  The  mean  minus  one  standard  deviation 
line  underestimates  the  area  underlain  by  coal  but 
approximates  it  much  more  closely  than  the  "best"  estimate. 


Experimental  variograms  computed  for  the  elevation  of 
the  Boundary  zone  (Figure  29)  are  similar  to  those  computed 
for  the  overlying  Estevan  zone  (Figure  24).  The  presence  of 
a  drift  is  indicated  by  the  divergence  of  the  directional 
variograms  beyond  a  distance  of  3  km.  The  variogram 
calculated  along  the  approximate  strike  (65°)  has  the  lowest 
drift  component  and  is  interpreted  as  the  variogram  of  the 
residual  component.  The  model  fitted  is  composed  of  nested 
spherical  structures  with  parameters! 

Ci  =  300  ft2  ai  =  500  metres 

C2  =  1200  ft2  a 2  =  15000  metres 

Although  a  drift  of  order  2  minimized  the  sum  of  squares  of 
the  residuals  (Table  2)»  the  result  produced  by  the  first 
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FIGURE  29.  BOUNDARY  ZONE  SURFACE,  VAR  I OGRAMS  AND  MODEL. 
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order  drift  is  similar*  The  first  order  drift  was  chosen 
primarily  because  it  permitted  extrapolation  over  a  wider 
area  but  secondarily  because  the  structural  analysis  for  the 
zone  may  have  been  biased  by  the  concentration  of 
observations  in  a  structurally  complex  area  near  subcrop* 


The  surface  of  the  B 
kriging  is  shown  in  Figur 
estimation  error  for  the 
in  Figure  30*  The  average 
between  15  feet  and  20 
points*  The  cluster  of  dr 
mile)  in  the  southeast  co 
the  standard  deviation  in 


oundary  zone  estimated  by  universal 
e  12*  The  standard  deviation  of  the 
surface  is  contoured  and  displayed 
value  of  the  standard  deviation  is 
feet  and  decreases  to  zero  at  data 
illing  at  400  metre  spacing  (1/4 
rner  of  Township  2y  Range  9  reduces 
the  area  to  less  than  15  feet* 


The  isopach  grid  of  the  interval  between  the  bedrock 
surface  and  the  Boundary  zone  (Figure  31)  was  created  by 
grid  manipulation*  The  gradient  on  the  isopach  surface  is 
greater  than  1°  near  the  predicted  subcrop  boundaries  except 
to  the  nor th— no r th wes t  under  the  city  of  Estevan  and  to  the 
southwest  of  the  bedrock  channel  in  Township  ly  Range  9*  The 
zone  does  not  approach  the  bedrock  surface  to  the  east  but 
is  bounded  by  a  lateral  change  to  sand*  The  standard 
deviation  of  the  estimation  error  for  the  isopach  is  shown 
in  Figure  32*  The  grid  values  average  between  20  to  25  feet 
of  standard  deviation  and  a  large  area  of  over  30  feet 
occurs  under  the  city  of  Estevan  in  a  data  void*  This  grid 


was 


manipulated  with  the  isopach  grid  (Figure  31)  resulting 
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FIGURE  30-  STANDARD  DEVIATION  OF  THE 

ESTIMATION  ERROR, 
BOUNDARY  ZONE  SURFACE. 
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in  the  predicted  location  of  Boundary  zone  subcrop  at  four 
probability  levels  (Figure  33)* 

The  "best"  estimate  of  subcrop  location  closely 
approximates  the  observed  borehole  results  except  southwest 
of  the  bedrock  channel*  Here»  geological  interpretation  of 
the  bedrock  topography  suggests  that  pre- glacial  erosion  has 
removed  the  Boundary  zone  along  a  nor th— northeast- trending 
valley*  The  Boundary  zone  intersected  in  boreholes  to  the 
northwest  of  this  valley  exists  probably  as  a  remnant 
isolated  by  erosion*  Detailed  drilling  at  400  metre  spacing 
in  the  southeast  quarter  of  Township  2t  Range  9  indicates 
that  the  Boundary  zone  is  structurally  disturbed  near  the 
subcrop  edge,  varying  more  than  30  feet  in  elevation  over  a 
distance  of  less  than  one  kilometer*  If  the  zone  is  net 
faulted,  dips  of  up  to  150  feet  per  mile  are  indicated  over 
a  distance  of  400  metres  and  nor thwes t— trend ing  swales  with 
an  amplitude  of  25  feet  spaced  at  about  500  to  800  metre 
intervals  may  be  interpreted*  Although  there  is  evidence 
that  the  Tertiary  coal  measures  are  faulted  in  the  area 


( Holter, 

1972, 

P* 

176  )»  the  problem 

sha  How 

boreh 

oles 

precluded  a 

interpretation* 

The 

presence  of  the 

f aul ted 


s true  tural 


in  boreholes  corresponds  roughly  to  structural  lows  in  the 
Boundary  zone* 


Detailed  drilling  indicates  that  the  structure  of  the 
coal  zone  near  the  subcrop  edge  is  locally  complex  due  to 
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glaciotec tonic  deformation  and  cannot  be  sampled  adequately 
at  a  drill  spacing  of  more  than  400  metres  nor  can  it  be 
adequately  represented  by  estimates  at  grid  nodes  spaced  500 
metres  apart*  The  local  structural  information  required  to 
predict  the  location  of  subcrop  in  the  Estevan  coalfield  can 
be  obtained  only  by  detailed  drilling*  Observation  points 
spaced  more  widely  than  the  dimension  of  local  structures 
will  sample  more  than  one  structure  and  a  false  trend  would 
be  Inferred  the  data*  Inadequate  sampling  results  in 
regional  extrapolation  of  false  trends  into  unsampled 
regions • 

The  uncertainty  of  correlation  compounded  by  a  complex 
local  structure  extrapolated  into  data  voids  accounts  for 
the  overe s t ima ti on  of  the  area  underlain  by  the  Estevan 
zone*  A  similar  although  much  more  minor  problem  occurs  for 
the  Boundary  zone  southwest  of  the  bedrock  channel*  The 
local  structure  is  reflected  in  the  variogram  models  of  the 
residual  component  for  the  coal  zone  surfaces  which  have  a 
short  range  (500  and  750  metres)  and  a  relatively  high  sill 
(225  ft2  to  300  ft2)*  The  grids  of  the  standard  deviation  of 
the  estimation  error  show  the  rapid  increase  in  uncertainty 
with  distance  away  from  a  borehole*  This  uncertainty 
flattens  out  beyond  about  1000  metres*  The  coarseness  of  the 
estimation  grid  (500  m)  relative  to  the  range  of  the 
variogram  (500  m  to  700  ra)  results  in  inadequate 
representation  of  the  narrow  area  surrounding  each  borehole 
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in  which  the  geology  is  known  with  grea 
of  the  crude  grid  approxiaation ,  b 
intersect  the  zone  lie  outside  the 
confidence  cumulative  probability  zones 
Figures  28  and  33 •  A  finer  grid  at  lea 
of  the  local  structure  (range  of  the 
required  for  adequate  representation 
variance*  The  cost  of  estimating  a  la 


nodes  may  be*  as  in  this  case,  prohibit 


t  confidence*  Because 
oreholes  which  do 
margins  of  the  higher 
(•84  and  *981)  in 
st  one  half  the  range 
variogram  model)  is 
of  the  estimation 
rge  number  of  grid 
ive  • 


An  alternative  approach  to  the  problem  of  subcrop 
estimation  is  conditional  simulation  (Chiles,  1976,  pp*  70— 
71)*  A  conditional  simulation  is  a  simulation  of  a 
regionalized  variable  with  the  following  characteristics: 

—  it  has  the  same  covariance  as  the  phenomenon 

—  it  passes  through  the  data  points* 

There  is  an  infinite  number  of  such  simulations  with  the 
same  covariance  which  pass  through  the  data  values*  If  the 
simulations  were  averaged,  the  average  value  at  any  point 
would  be  the  value  estimated  by  kriging  while  the  variance 
would  be  the  estimation  variance*  Conditional  simulations 
may  be  more  useful  than  kriging  estimates  since  the 
simulation  values  show  the  true  variability  of  variable 
values  whereas  the  kriging  estimates  are  smoothed*  The 
conditional  simulations  may  be  considered  as  possible  models 
of  the  true  (unknown)  surface  and  the  one  which  seems  to  fit 


the  other  information  (such  as  boreholes  which  did  not 
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Intersect  a  zone)  can  be  selected*  If  the  cost  of  creating 
conditional  simulations  at  appropriate  scale  is  not  too 
great*  the  method  could  be  a  viable  alternative  to  the 
kriging  technique* 


4  • 2  £rp^abj>.\jty  sl£  JttlS  Area  Underlain  hx  Coaj. 

Most  damaging  of  all  criticisms  of  area  estimation  by 
subcrop  lines  is  that  lithologic  boundaries  may  account  for 
much  of  the  perimeter  of  an  area*  These  boundaries  are 
either  depositlonal  or  erosional  in  nature*  The  depositional 
extent  of  coal  is  controlled  by  the  paleoenvironment ;  a 
complex  combination  of  factors  including  climate* 
topography*  hydrology  and  time*  Reconstruction  of  the 
lateral  facies  relationship  of  a  peat*f ormi ng  environment  is 
difficult  from  the  limited  exposures  available*  McCammon 
(1975)  demonstrated  that  even  the  use  of  sophisticated 
sampling  rules  devised  with  prior  knowledge  of  the  facies 
distribution  cannot  provide  significant  improvements  over  a 
systematic  sample  pattern  in  reconstruction  of  the  lateral 
facies  relationships  in  the  Mississippi  delta*  A  sinuous 
erosional  channel  is  an  example  of  a  facies  which  is 
difficult  to  trace  in  subsurface* 


Because  the  prediction  of  lithogic  boundaries  requires 
a  great  deal  of  interpretation,  it  is  not  amenable  to  a 
statistical  approach  based  on  objective  probability*  The  use 
of  subjective  probability  and  cumulative  probability 
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distributions  of  variables  in  resource  potential  evaluation 
has  been  described  in  EMR  (  1977  )  and  Newendorp  (  1975  )•  The 
variable  under  study  here  is  the  area  underlain  by  coal* 
Three  parameters  of  Its  cumulative  probability  distribution 
were  estimated:  the  minimum  probable  area  underlain  by  coal 
(  P[  X  >  x]  =  1  )  *  the  maximum  probable  area  (  P[  X  >  x]  =  0)  and 
the  most  likely  area  CP[  X  =  x]  is  maximum)*  These  values 
were  estimated  by  outlines  of  the  maximum*  minimum  and  most 
likely  extent  of  the  coal  zones  under  the  most  optimistic* 
pessimistic  and  reasonable  interpretations  of  all  the 
available  data*  These  estimates  are  shown  for  the  Estevan 
and  Boundary  zones  in  Figure  34  and  Figure  35  and  the  values 
are  tabulated  in  Table  3*  These  values  were  interpolated  to 
form  a  continuous  cumulative  probability  distribution 
(Figure  36)  using  the  triangular  distribution* 


la  Els  2 


Area  estimates  of  the  coal  zones*  Estevan  coalfield 
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Figure  36  -  cumulative  probability  distributions ,  areal  extent  of  the  coal  zones 
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4 • 3  Aygfftge  Xhiskasss  12JC  Net  Coal  jjj  a  ^one 

The  average  thickness  of  net  coal  was  estimated  by 
kriging  within  the  best  estimate  (modal)  of  the  areal  extent 
of  the  zones  (Figures  34  and  35 )« 

The  experinental  variograms  and  sample  variance  for  the 
net  coal  in  the  Estevan  zone  are  shown  in  Figure  37*  The 
directional  variograms  do  not  differ  markedly  and 
directional  anisotropy  is  not  indicated*  The  model  fitted  is 
a  combination  of  two  isotropic  spherical  schemes  with  a 
nugget  effect*  The  parameters  of  the  model  are: 

C0  =  1.0  ft2 

C 1  =  4*0  ft2  a i  =  1000  metres 

C 2  —  5*5  ft2  a 2  —  4000  metres 

The  nugget  effect  is  used  to  represent  a  very  short  range 
component  of  the  uncertainty  caused  by  the  inclusion  or 
rejection  of  thin  overlying  or  underlying  seams  into  the 
zone  according  to  the  selection  criteria*  The  variogram 
model  rises  quite  rapidly  within  the  first  kilometre  and 
then  levels  out  at  about  4  kilometres*  The  average  variogram 
remains  at  this  sill  until  drift  effects  cause  a  rise 
slightly  beyond  9  kilometres*  The  drift  effect  is  not 
important  and  was  not  modelled* 

The  variance  of  the  area  and  the  average  covariance  of 
the  samples  to  the  area  were  estimated  as  discrete  sums  over 
a  500  metre  grid  of  1789  points  representing  an  area  of 
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DISTANCE  .  METRES 

FIGURE  37.  ESTEVRN  NET  COAL  THICKNESS.  VRRIOGRAMS  AND  MODEL. 
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about  170  square  miles*  Sample  values  fron  160  boreholes 
formed  the  basis  of  the  average  thickness  estimate*  Since 
the  Estevan  zone  variogram  model  is  composed  of  a  short 
range  component  <1000  »)  accounting  for  about  40%  of  the 
variation,  the  choice  cf  500  metres  as  a  grid  size  resulted 
in  about  a  5%  overestimation  of  the  estimation  variance*  As 
grid  size  was  decreased,  the  weights  calculated  by  the 
kriging  system  did  not  change  markedly  nor  did  the  average 
thickness  but  the  variance  of  the  area  decreased  towards  a 
stable  value* 


The  average  net  coal  thickness  of  the  Estevan  zone  is 
estimated  to  be  8*9  feet  and  the  estimation  variance  is 
0*085  ft2*  This  compares  to  an  unweighted  mean  of  9*5  feet 
and  variance  of  the  mean  of  0*077  ft2*  The  unweighted  mean 
appears  biased  on  the  high  side  by  about  7* 


The  experimental  variograms  and  va 
coal  thickness  of  the  Boundary  zone  is 
The  net  coal  thickness  appears  to 
nugget  effect*  A  transition  phenomenon 
(1000  m)  was  modelled  within  the  najcr 
The  parameters  of  the  nested  spherical 


Cl 

=  1.0 

ft2 

ax  =  1000 

metre  s 

c2 

=  4.5 

ft2 

a  2  =  6000 

me  t  re  s 

The  average  net  coal  thickness 


riance  for  the  net 
plotted  in  Figure  38* 
be  isotropic  with  no 
with  a  short  range 
long  range  structure* 
models  are* 


over  the  area  was 


estimated  from  161  sample  values  over  an 
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DISTANCE  .  METRES  X104 

FIGURE  38.  BOUNDARY  NET  COAL  THICKNESS.  VAR  I QGRAMS  AND  MODEL. 
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square  miles  represented  by  2115  points  located  on  a  500 
metre  grid*  Because  of  the  discrete  summation  of  the  area 
variance  and  average  covariances*  the  estimation  variance 
may  be  3%  to  4%  too  high*  The  relative  error  is  less  for  the 
Boundary  zone  than  for  the  Estevan  zone  because  the  short 
range  structure  is  a  less  significant  component  of  the 
Boundary  zone  model* 

The  average  thickness  estimate  is  3*8  feet  and  the 
estimation  variance  is  0*057  ft2*  The  unweighted  average  of 
sample  values  is  5*3  feet  and  the  variance  of  the  mean  is 
0*043  ft2*  The  unweighted  average  is  biased  towards  a 
greater  thickness  because  of  the  disproportionate  number  of 
samples  taken  near  the  subcrop  line  to  the  northwest  (Figure 
11)*  In  the  case  of  the  Boundary  zone*  the  unweighted 
estimate  is  over  a  third  too  highf  demonstrating  the 
importance  that  spatial  correlation  has  in  the  selection  of 
an  appropriate  unbiased  estimator  for  the  net  coal 
thickness • 

4.4  langagg  SalitBales 

Having  estimated  and  represented  the  component 
variables  of  the  volume  equation  (average  thickness  and 
area)  in  the  form  of  cumulative  probability  distributions* 
the  addition  of  a  tonnage  tonnage  factor  based  upon  specific 
gravity  permitted  the  calculation  of  in— situ  coal  tonnages* 


The  mean  specific  gravity  was  estimated  to  be  1*24  with  a 
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minimum  of  1*19  and  a  maximum  of  1*29  (  J*  Irvine*  personal 
communication)*  These  values  were  taken  as  ±2  limits  of  a 
normal  distribution*  The  cumulative  normal  distribution  was 
used  for  both  the  tonnage  factor  and  the  average  thickness* 

The  evaluation  of  the  resource  equation  was  carried  out 
by  a  simulation  technique  described  by  Roy  (1975)*  In  this 
process  component  variable  distributions  are  sampled 
independently  and  the  sample  values  are  multiplied  to  form 
one  estimate  of  the  resource  tonnage*  The  simulated  sampling 
process  was  repeated  5000  times  and  the  results  were 
accumulated  as  a  coal  tonnage  cumulative  probability 
distribution*  The  simulation  technique  is  particularly 
useful  if  the  distributions  of  the  variables  in  the  resource 
potential  equation  are  irregular  or  highly  skewed  and  cannot 
be  represented  by  analytical  expressions*  The  area  underlain 
by  coal  is  an  example  of  an  assymetric  distribution* 


The  tonnage  distributions  created  by  the  simulation 
runs  for  the  Estevan  and  Boundary  zones  are  displayed 
graphically  in  Figure  39*  The  parameters  of  the  tonnage 
distributions  and  those  of  the  component  variables  are 
tabulated  in  Table  4*  The  resource  quantity  estimated  for 
the  Estevan  zone  is  twice  that  of  the  Boundary  zone*  The 
coefficient  of  variation  and  max lroum/mean  ratio  of  tonnages 
suggest  that  the  Estevan  zone  is  estimated  with  narrower 
relative  confidence  limits  than  the  Boundary  zone* 
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Parameters  of  estimated  coal  resource  quantity 
distributions  and  component  variables 


Distr ibuti on 

( Units  ) 

Mean 

Vari ance 

Coefficient 

o  f 

variation 

maximum 

mean 

EstSYftQ 

tonnage 
<  tons  ) 

l  • 59*1 09 

•200*10 1 7 

8.9 

1.3 

area 

( sq •  mi les ) 

165.5 

175.7 

8.0 

1.2 

thickness 
(  feet  ) 

8.92 

•  0854 

3.3 

1.08 

tonnage 

factor 

( tons/mi2f t  ) 

1 • 08*106 

•473*106 

2.0 

1.04 

Boundary 

tonnage 
( tons  ) 

796* 1 06 

•  812*10*  6 

11.3 

1.4 

area 

(  sq •  miles  ) 

192.4 

309.7 

9.1 

1.2 

thickness 
( fee  t  ) 

3.83 

.05735 

6.25 

1.14 

tonnage 

factor 

(  tons/mi2f t  ) 

1 • 08*106 

•473*106 

2.0 

1.04 

Egteyflu 

and  Boundary 

tonnage 
(  tons  ) 


2.40*1Q9 


278*10*  7 


6.9 


1  .3 
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If  the  confidenc 


are  taken  to  be  .977  th 


the  mean  ±2  standard 


mean  ±18%  for  the  Estev 


Boundary  zone.  The  Est 


fall  in  the  measured  re 


while  the  Boundary  r 


outside  the  ±20%  cut-of 


e  limits  of  the  tonnage  distribution 
en  the  range  of  values  considered  is 
deviations.  This  corresponds  to  the 
an  zone  and  the  mean  ±23%  for  the 
evan  resource  quantity  estimate  would 
source  category  of  McKelvey  (  1976  ) 
esource  quantity  estimate  falls  just 
f  • 


The  percentage  contribution  of  the  component  variables 
to  the  variance  of  the  tonnage  distributions  is  given  in 
Table  5.  In  both  zones  the  uncertainty  in  the  areal  extent 
of  the  zone  accounts  for  over  half  the  total  uncertainty. 
The  contribution  of  the  average  thickness  variable  ranks 
second  and  is  more  significant  for  the  Boundary  zone.  The 
coefficient  of  variation  for  thickness  in  the  Boundary  zone 
(Table  4)  is  almost  twice  that  in  the  Estevan  zone.  The 
coefficient  of  variation  for  the  area  of  the  Boundary  zone 
is  only  slightly  greater  than  that  of  the  Estevan  zone.  In 
neither  case  does  the  variance  of  the  tonnage  factor 
contribute  more  than  6%. 


The  thickness  of  the  Estevan  zone  is  apparently  well 
estimated  from  the  fairly  regular  sample  distribution 
(Figure  5).  The  area  estimate  suffers  from  the  long  subcrop 
perimeter  which  are  structurally  complex  and  inadequately 
sampled  and  from  the  presence  of  an  erosional  channel,  the 


location  and  extent  of  which  is  difficult  to  define  (Figure 
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Components  of 

the  variance,  tonnage 

distributions 

Z  one 

Es  tevan 

Boundary 

area 

SI.  1 

65.8 

thickness 

13.8 

31.0 

tonnage  factor 

5.  1 

3.2 

tonnage 

100*0  percent 

100.0  percent 
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34)*  The  thickness  of  the  Boundary  zone  is  less  well  defined 
because  of  the  non— sy s temat ic  location  of  samples  (Figure 
11  )  and  because  the  eastern  boundary  of  the  areal  extent  is 
subject  to  several  interpretations  which  vary  widely  (Figure 
35  )•  In  both  cases*  any  significant  reduction  in  the 
uncertainty  is  dependent  upon  more  detailed  drilling  to 
define  the  boundaries  of  the  coal  zone*  Additional  boreholes 
in  data  voids  would  improve  the  precision  of  the  thickness 
estimate*  particularly  for  the  Boundary  zone  in  which  more 
than  80%  of  the  variability  occurs  within  a  range  of  one  to 
six  kilometres  of  a  point*  The  location  of  one  additional 
borehole  which  would  maximize  the  reduction  in  the 
estimation  variance  of  the  average  thickness  could  be 
determined  using  methods  described  in  Delfiner  and  Delhomme 
(  1975  ). 
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CHAPTER  5 


5. 


Ajaa.ilsallgns  la  S.ygt emails  Sagaling 


■Small  Ss.al^  Area 


Three  sections  in  the  small  scale  area  (Figure  18)  were 
drilled  over  a  period  of  three  years  at  grid  spacings  of 
800'  400  and  200  metres  with  some  infill  drilling  in  the 
fourth  year  at  100  metre  spacing  to  define  an  erosional  area 
(Figure  14)*  The  proportion  of  boreholes  which  intersect 
coal  is  an  unbiased  estimator  of  the  proportion  of  the  area 
underlain  by  coal  and  the  arithmetic  mean  of  the  net  coal 
thickness  in  holes  penetrating  the  zone  is  an  unbiased 
estimator  of  the  average  net  coal  thickness* 


The  precision  of  a 
coal  is  a  function  o 
thickness*  the  shape  of 
between  observations* 
estimate  adopted  in  thi 
standard  deviations  exp 
of  the  estimate: 


n  estimate  of  the  volume  of  in— situ 
f  the  natural  variability  of  the  zone 
the  coal  area  and  the  distance 
The  measure  of  precision  of  an 
s  study  is  defined  to  be  two  relative 
ressed  as  a  percentage  of  the  value 


P  =  100*(  2*S)/m  (  5.1  ) 
where  S2  is  the  estimation  variance  and  m  is  the  value  of 
the  estimate*  If  the  error  made  in  estimating  the  true 
unknown  M  by  m  is  normally  distributed  M  occurs  within  the 
interval  m±2S  with  confidence  0*98* 
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The  relationship  between  the  precision  of  estimation 
and  the  amount  of  exploration  effort  is  not  a  simple  one* 
Once  knownt  the  relationship  can  be  used  to  select  the  drill 
spacing  appropriate  to  an  acceptable  precision  level,  and 
the  rate  of  improvement  in  precision  as  a  function  of 
exploration  effort  can  also  be  studied* 

A  simplistic  model  of  the  relationship  between  the 
precision  of  the  estimate  (P)  and  the  number  of  boreholes 
(N)  drilled  in  a  limited  area  is 

InP  =  c-klnN  (5*2) 

where  c  and  k  are  constants*  If  k  equals  one,  the  model  may 
be  rewritten  simply  as 
P  =  c/N 

and  doubling  the  exploration  effort  effectively  halves  the 
uncertainty*  If  k  is  greater  than  one,  the  return  on  the 
increased  effort  is  magnified  whereas  if  k  is  less  than  one, 
the  increase  in  precision  follows  a  law  of  diminishing 
returns*  In  practice,  k  may  also  be  a  function  of  N, 
increasing  as  the  information  provided  by  more  samples  has  a 
synergistic  effect*  The  value  and  behaviour  of  k  can  be 
examined  by  plotting  InP  against  InN  for  experimental  data* 
The  slope  or  tangent  to  the  curve  Joining  the  points  is  the 


value  of  k* 
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5*1  £r.egislpn  £>1  Area  fslircale 


On  the  basis  ef  a 
study  area  is  composed 
1*  an  area  in  which  th 
2*  an  area  in  which 
eroded  at  the  bedro 
3 •  an  area  in  which  th 
deposi ted • 

The  narrow,  irregular 
problem  in  making  the  d 
incomplete  sections  o 
metres  or  more  apart# 
precision  at  all 
intersecting  any  coal  i 
intersection  and  the 
coal  and  non— coal* 


11  the  available  information,  the 
of  three  areas  (Figure  18): 
e  complete  coal  zone  exists* 
the  upper  seams  of  the  zone  have  been 
ck  unconformity* 

e  zone  has  been  eroded  or  was  not 

width  of  the  erosional  edge,  posed  a 
istinction  between  the  complete  and 
f  coal  zone  in  boreholes  spaced  400 
For  the  purpose  of  comparing  the 
observation  distances,  a  borehole 
n  the  zone  was  considered  a  coal 
area  was  divided  into  only  two  areas: 


The  area  of  non— coal  was  calculated  as  the  product  of 
the  number  of  boreholes  which  did  not  intersect  coal  and  the 
area  of  a  square  with  sides  the  length  of  one  observation 
spacing*  Boreholes  along  the  edge  and  at  the  corners  of  the 
study  area  contributed  0*6  and  0*25  respectively  of  the  area 
of  an  interior  borehole*  The  relative  variance  of  the  non- 
coal  area  estimate  was  calculated  using  Natheron's 
approximation  formula  (3*20  in  chapter  3)*  The  estimation 
variance  and  precision  for  the  area  underlain  by  coal  was 


derived  from  the  relative  variance  of  the  non— coal  area 
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The  parameters  of 
and  estimation  varianc 
listed  in  Table  6.  The 
the  number  used  in  a  c 
fact  achieved  at  the  20 
the  number  of  cells  i 
interpreted  from  Figure 
because  of  the  arbitra 
Nt  and  N2  are  one— half 
in  the  directions  of  th 
fractions  because  of  th 


the  formulaf  the 


e  values  for  thr 


number  of  boreh 


omplete  grid  pat 


0  metre  spacing. 


n  which  coal  was 


18*  It  contains 


ry  boundaries  of 


the  estimated  le 


e  grid  lines 


e  corner  and  edg 


resulting  precision 
ee  drill  spacings  arc 
oles  corresponds  to 
tern  which  was  not  in 
The  parameter  N  is 
not  observed  and  was 
a  decimal  fraction 
the  area*  Parameters 
ngth  of  the  perimeter 
and  contain  decimal 
e  boreholes* 


Matheron  (1971*  pp*  27—28)  cautions  that  formula  3*20 
expresses  the  estimation  variance  only  if  the  spacing  of 
observations  is  small  and  that  it  neglects  a  fluctuating 
term  which  may  be  quite  significant*  This  fluctuating  term 
results  from  the  unique  combination  of  the  random  location 
of  the  grid  origin  and  the  frequency  of  observations 
relative  to  the  frequency  of  boundaries  between  coal  and 
non— coal  areas*  Referring  to  Figure  18*  it  appears  that  an 
800  metre  grid  spacing  is  larger  than  the  average  diameter 
of  the  non— coal  areas*  Because  an  assumption  basic  to  the 
calculation  is  violated,  the  precision  of  the  800  metre 


estimate  should  be  taken  as  an  approximation* 


By  way  of  comparison* 


the  precision  of  the  area 
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Precision  of  area  estimation  using  Matheron* s  approximation 
formula  for  relative  estimation  variance* 


Drill 
spacing 
( metres  ) 

800 

400 

200 

Number  of 
holes 

21 

65 

225 

Ni 

1.75 

7 

17 

n2 

1.75 

7 

14 

N 

1.75 

9.25 

32.75 

area  of  coai. 

( sq •  miles) 

2.56 

2.42 

2.48 

variance 
(  ft2  ) 

.  0256 

.0064 

.0009 

precision 

25% 

6.5% 

2.4% 

precision 
( binomial  ) 

19% 

11% 

5.9% 
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estiaate  may  also  be  calculated  based  on  an  assumption  that 
areal  proportion  estimates  are  distributed  as  the  binomial 
distribution*  These  precision  estimates  are  included  in 
Table  6*  The  precision  data  in  Table  6  have  been  plotted 
against  the  number  of  observations  on  full  logarithmic  axes 
in  Figure  40* 

According  to  Matheron's  approximation  formula,  the 
uncertainty  decreases  more  rapidly  than  the  grid  spacing  but 
that  it  does  not  decrease  as  rapidly  as  the  number  of 
observations  increases*  The  slope  of  the  line  segment 
between  200  n  and  400  m  is  about  -'•>0*8  which  is  an 
improvement  over  the  constant  slope  of  —0*5  for  the  binomial 
model*  The  slope  of  the  curve  for  the  approximation  formula 
can  be  expected  to  be  close  to  unity  but  will  vary  since  the 
relationship  also  depends  on  the  parameters  Nj,  and  N2  which 
express  the  complexity  of  the  coal  zone  boundary  apparent  at 
various  sample  spacings* 


5.2  structural  Analysis  Eaplaraliflii  Jl^ta 


Precision  estimates  for  the  average  net  coal 
are  based  on  the  model  f i tted  to  experimental 
The  accuracy  of  the  precision  estimates  depends 
recovery  of  the  underlying  variogram  model  from 
data*  If  the  raw  sample  data  contain  more 
population,  the  natural  variability  of  the 


thickness 
variograms • 
upon  the 
the  sample 
then  one 
net  coal 


thickness  may  be  obscured  by  the  effect  anomalous  outlying 
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FIGURE  40- RELATIVE  PRECISION  OF  THE  AREA 
ESTIMATE  AS  A  FUNCTION  OF  THE  NUMBER  OF 
OBSERVATIONS.  SMALL  SCALE  AREA. 
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values  have  on  the  experimental  variogram*  That  small 
samples  often  do  not  provide  stable  estimates  of  the 
population  meant  variance  and  variogram  is  borne  out  by 
Table  7  in  which  the  sample  mean  variance  and  precision  for 
net  coal  thickness  are  calculated  for  three  sample  spacings 
using  the  available  data  in  the  study  area*  The  calculations 
were  made  for  all  boreholes  together  and  also  for  boreholes 
which  show  only  a  complete  zone  section*  The  number  of 
observations  in  the  truncated  areas  is  also  tabulated*  The 
mean  appears  to  be  a  relatively  stable  estimate  whereas  the 
variance  decreases  as  the  number  of  observations  increases* 
The  variance  increases  for  all  boreholes  at  all  spacings 
(Table  7*  column  4)  because  of  the  effect  of  additional 
holes  penetrating  incomplete  sections  of  the  coal  zone* 

As  a  homogeneous  population*  the  variance  for  the 
complete  zone  is  smaller  than  that  for  the  mixed  populations 
at  all  sample  spacings*  The  precision  of  estimated  average 
net  coal  thickness  of  the  complete  zone  is  about  75%  of  the 
mixed  population  and  the  mean  is  5%  higher*  Table  8 
demonstrates  that  similar  results  hold  for  all  the  samples 
collected  in  the  small  scale  area*  With  an  increased  number 
of  samples*  the  precision  values  for  the  entire  small  scale 
area  are  lower  than  for  the  study  area  although  the  mean  and 
variance  are  comparable* 

The  sample  variance  perceived  at  low  information  levels 
is  higher  than  the  real  population  variance  and  leads  to 
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Sample  mean  and  variance  of  net  coal  thickness  observations 
at  three  drill  spaclngs.  Data  are  for  three  sections  only# 


Drill  spacing* 

( metres  ) 

800 

400 

200 

All 

*  Complete  Zone  * 

Number  of 
observations 

18 

44 

157 

173 

Mean  (  f t ) 

12.7 

12.6 

12.6 

12.6 

Variance  (  ft2 ) 

9.3 

5.0 

4.25 

4. 16 

Precisi on  ( % ) 

11.3 

5.35 

2.60 

2.46 

*  Incomplete  Zone 

* 

Number  of 
observations 

3 

7 

22 

4  t 

*  All 

Boreholes  * 

Number  of 
observations 

21 

51 

179 

217 

Mean  (  f t ) 

12. 1 

12.0 

12.0 

11.65 

Variance  (  ft2  ) 

14.9 

8.4 

7.47 

8.55 

Precision  ( %  ) 

13.9 

6.8 

3.4 

3.4 
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Sample  mean  and  variance  of  net  coal  thickness  observations 
at  three  spacings.  All  observations  in  the  small  scale  area* 


Drill  spacing 
( metres  ) 

800 

400 

200 

All 

*  Complete  Zone  * 

Number  of 

observations 

62 

168 

344 

360 

Mean  ( f t  ) 

11.6 

12.0 

12.2 

12.2 

Variance  ( ft2  ) 

7.9 

4.76 

4*23 

4. 18 

Preci si  on  (%) 

6.1 

2.8 

1.8 

1.8 

*  Incomplete  Zone 

$ 

Number  of 

observations 

7 

12 

28 

50 

Mean  (  f t  ) 

7.47 

Variance  (  ft2  ) 

9.7 

Preci si on  ( %  ) 

* 

All  Coal  * 

11.3 

Number  of 
ob  s  e  rv  a  t  i  o  n  s 

69 

180 

372 

410 

Mean  (  f t  ) 

12.0 

11.6 

i  1.78 

11.63 

Variance  (ft2) 

12.3 

7.4 

6.64 

7.25 

Preci si on  ( %  ) 

7.0 

3.5 

2.3 

2.3 

' 
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precision  estimates  which  are  higher  than  those  based  on  the 
population  variance*  Precision  estimates  made  with  only  a 
few  samples  provide  conservative  confidence  limits  for 
resource  estimates* 


Experimental  variograms  were  constructed  using  data 
from  boreholes  selected  with  complete  sections  of  the  coal 
zone  and  using  data  from  all  boreholes*  The  data  available 
at  the  end  of  each  year  of  drilling  were  used  to  illustrate 
the  effect  of  the  grid  spacing  on  perception  of  the 
underlying  variogram*  Variograms  of  the  selected  samples 
represent  the  natural  variability  of  the  coal  zone  whereas 
variograms  for  all  samples  contain  a  large  variance 
component  for  the  small  number  of  samples  from  areas  where 
the  coal  was  truncated  by  erosion* 


Variograms  for  the  first  year  (Figures  41  and  42)  do 
not  show  any  significant  spatial  structure  and  would  be 
interpreted  as  the  random  type  (Figure  21)*  Variograms 
calculated  after  the  second  year  (Figure  43  and  44)  indicate 
some  spatial  structure  and  a  level  of  variance  considerably 
reduced  from  the  previous  year*  The  variogram  of  the 
selected  samples  (Figure  43)  has  a  possible  short  range 
structure  which  rises  to  a  sill  of  3*75  ft2  at  about  600  to 
800  metres*  The  directional  variogram  for  an  azimuth  of  45° 
rises  above  the  sample  variance  beyond  2500  metres  in 
response  to  a  thinning  of  net  coal  to  the  southwest  (Figure 
19).  Variograms  based  on  all  the  samples  (Figure  44) 
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DISTANCE  ,  METRES 

FIGURE  41.  AVERAGE  VARIOGRAM  OF  NET  COAL  THICKNESS,  SELECTED  SAMPLES  -  YEAR 
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0ISTRNCE  ,  METRES 

FIGURE  42.  AVERAGE  VAR  I OGRAM  OF  NET  COAL  THICKNESS,  ALL  SAMPLES  -  YEAR 
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DISTANCE  ,  METRES 

FIGURE  43.  VflR I OGRAMS  OF  NET  CGRL  THICKNESS.  SELECTED  SAMPLES  -  YEAR 
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DISTANCE  ,  METRES 

FIGURE  44.  VRRIOGRAMS  OF  NET  COAL  THICKNESS.  ALL  SAMPLES  -  YEAR 
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demonstrate  the  effect  of  outlying  values  from  the  partly 
eroded  coal  zone*  The  directional  variograms  show  a  greater 
scatter  around  the  average  variogram  and  a  linear  increase 
in  the  variogram  with  distance* 

Variograms  calculated  from  selected  samples  after  the 
third  year  of  drilling  (Figure  45)  are  almost  flat  and  the 
drift  in  the  45°  direction  is  defined  more  clearly* 


Additional  samples 

drilled 

at 

100 

metre  spacing  in 

the 

fourth  year  add  very 

little 

to 

the 

variogram  model 

for 

selected  samples  (Figure  46 )«  A  structure  with  a  range  of 
about  800  metres  is  suggested  in  both* 

When  all  samples  are  considered  in  the  variogram 
calculation  after  the  third  year,  the  shape  of  the  resulting 
variograms  ( Figure  47  )  parallels  those  for  the  selected 
samples  (Figure  45)  at  a  level  of  variance  about  2  ft2 
higher*  A  short  range  structure  with  a  range  of  about  800 
metres  appears  to  be  the  most  likely  interpretation*  As 
expected,  the  variogram  of  samples  from  the  erosional  edge 
(Figure  48)  shows  no  spatial  structure* 

These  results  demonstrate  that  an  experimental 
variogram  calculated  from  a  small  number  of  samples  taken  at 
widely— spaced  locations  rarely  indicates  the  underlying 
structure  of  the  net  coal  thickness*  A  stable  estimate  was 
obtained  with  no  fewer  than  50  observations  located  with  a 
minimum  spacing  of  one— half  the  range  of  the  smallest 
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FIGURE  46.  VARIOGRAMS  OF  NET  CORL  THICKNESS*  SELECTED  SAMPLES  YEAR 
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FIGURE  47.  VARIGGRAMS  OF  NET  COAL  THICKNESS.  ALL  SAMPLES  -  YEAR 
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FIGURE  48.  AVERAGE  VAR  I OGRAM  OF  NET  COAL  THICKNESS.  SAMPLES  FROM  ERODED  AREAS 
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significant  structure*  If  observations  from  two  populations 
are  intermixed)  the  experimental  variogram  will  overestimate 
the  uncertainty  in  the  individual  populations  and  the  range 
of  the  major  phenomenon  will  be  obscured* 


Techniques  which  optimize  the  spacing  of  observations 
for  a  desired  precision  of  the  estimate  are  based  upon 
foreknowledge  of  the  underlying  variogram*  Unless  this 
variogram  can  be  identified  from  outcrop  exposures  or  mining 
data)  such  techniques  are  of  no  value  in  the  reconnaissance 
stage  of  regional  coal  resource  evaluation  programs*  The 
requirement  in  number  and  spacing  of  observations  for  a 
stable  variogram  estimate  may  not  be  met  until  the  pre- 
production  stage  of  exploration*  This  limitation  could  be 
overcome  if  the  known  variogram  of  an  analogous  well- 
explored  coal  zone  is  used  as  a  model  for  the  relatively 
unexplored  coal  zone* 

The  smallest  separation  of  sample  pairs  which  provides 
a  stable  variogram  estimate  is  about  200  metres*  Structures 
in  the  net  coal  thickness  which  have  an  average  diameter 
less  than  200  metres  cannot  be  identified  on  the  variograms* 
Because  the  behaviour  of  the  variogram  near  the  origin  has  a 
significant  effect  on  the  calculated  estimation  variance, 
some  assumptions  were  made  about  the  behaviour  of  net  coal 
thickness  in  very  local  areas*  The  most  reasonable 
assumption  geologically  Is  that  net  coal  thickness  is 


continuous  and  that  there  exists  a  structure  (or  structures) 


. 
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with  a  range! s)  shorter  than  200  metres*  The  implications  of 
the  choice  of  a  range  for  the  short  scale  structure  to  the 
precision  of  estimates  will  be  discussed  in  the  next 
section • 


The  interpretation  of  net  coal  thickness  as  a 
continuous  phenomenon  is  illustrated  for  the  average 
variograu  of  all  samples  after  the  third  year  of  drilling 
(Figure  49)*  The  parameters  of  this  model  are. 

Ci  =  3  ft2  ai  =  200  metres 

Cj  =  3  ft2  aj  =  800  metres 

The  choice  of  the  model  for  all  samples  after  the  third  year 
for  subsequent  precision  estimates  was  guided  by  two 
considerations: 

1*  The  precision  estimates  for  the  average  net  coal 
thickness  must  be  directly  comparable  to  results  derived 
for  the  area  estimates*  Thus  all  samples  representative 
of  the  coal  zone  had  to  be  included* 

2*  This  model  was  more  efficient  computationally  than  the 
alternative  of  modelling  the  selected  coal  zone  samples 
and  samples  from  the  erosional  areas  separately  and 
combining  the  variances* 


5.3  Precision  Q.1  t he  Ihigjingsa  EglljULl&S 

The  precision  of  the  average  net  coal  thickness  was 
calculated  for  nine  sample  grids  using  two  different 
assumptions*  The  first  is  based  on  the  assumption  that 
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observations  are  statistically  independent  and  the 
population  has  a  mean  of  12  feet  and  a  variance  of  6  ft2. 
The  precision  in  this  case  is  an  inverse  function  of  the 


square  root  of  the  number  of  samples: 
P  =  100  *  2  * 

1  N 


(5.3  ) 


The  second  type  of  precision  calculation  was  based  on  the 
assumption  that  the  net  coal  thickness  has.  a  variogram 
composed  of  a  two  nested  structures?  a  short  range  structure 
with  a  sill  of  3  ft2  at  200  metres  or  less  and  a  long  range 
structure  with  a  range  of  800  metres  and  a  sill  of  3  ft2. 
The  variance  of  this  model  is  the  sum  of  the  sill^  of  the 
component  models  6  ft2  and  the  population  mean  was  assumed 
to  be  12  feet.  The  variance  ( S2  )  of  the  error  made  in 
extending  the  value  of  a  central  borehole  to  a  square  of 
influence  the  size  of  the  grid  spacing  was  calculated  using 
auxiliary  functions  (Clark,  1976).  If  the  errors  made  in 
each  grid  cell  are  independent,  the  precision  of  the 
estimate  is  calculated  as  follows: 


P  =  100  *  2  *  _/  S2 

N 


(  5.4  ) 


The  extension  variances  ( S2 )  were  calculated  for  five 
ranges  assumed  for  the  short  range  structure  are  tabulated 
in  Table  9.  Varying  the  range  from  zero  to  200  metres  has 
the  effect  of  reducing  the  estimation  variance.  The  effect 
is  most  pronounced  at  grid  spacings  shorter  than  400  metres 
illustrating  that  the  choice  of  a  model  near  the  origin  is 
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Iflble  2. 


Variance  ( It2  )  of  the  error  made  In  extending  the  value  of  a 
central  borehole  to  a  square  L  metres  by  L  metres-* 
Calculations  were  performed  for  various  assumed  values  for 
the  range  of  the  variogram  short  range  structure* 


Range  of 
Variogram  (m) 

0 

10 

50 

100 

200 

i,  C  m  ) 

50 

3.  069 

2.985 

1.290 

0.633 

0.345 

100 

3.  137 

3.118 

2.549 

1.359 

0.7  02 

200 

3.276 

3.271 

3.143 

2.688 

1.498 

300 

3.418 

3.416 

3.361 

3. 173 

2.356 

400 

3.  564 

3.563 

3.533 

3.  431 

2.976 

600 

3.876 

3.875 

3.862 

3.819 

3.631 

800 

4.222 

4.221 

4.214 

4.190 

4.  089 

1200 

4.938 

4.938 

4.935 

4.925 

4.881 

1600 

5.412 

5.412 

5.410 

5.404 

5.380 

« 
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significant  in  the  calculation  of  precision  and  particularly 
the  improvement  of  precision  with  increased  exploration 
effort#  The  precision  of  the  average  thickness  estimate  was 
calculated  for  three  error  variance  models:  statistical 
independence*  the  0  metre  range  and  the  200  metre  range 
variogram  models  (Table  10  )•  Precision  values  for  variogram 
models  with  ranges  between  0  m  and  200  m  fell  between  the 


values  tabulated  for  0  and  200  metres# 


The  logarithm  of  precision  was  plotted  against  the 
logarithm  of  the  number  of  observations  from  Table  10  in 
Figure  50#  The  relationship  between  precision  and  the  number 
of  observations  in  the  case  of  statistical  independence  is 
strictly  linear  with  a  slope  of  -0#5#  The  variogram  with  the 
nugget  effect  represents  improved  precision  but  also 
displays  a  linear  relationship  with  slope  ■“0#5#  The  most 
continuous  variogram  of  net  coal  thickness  (200  metre  range) 
results  in  precision  which  is  not  a  linear  function  of  N  but 
improves  at  grid  spacings  of  400  metres  and  less#  The  value 
of  k  at  200  metres  is  about  -0 • 8  indicating  a  substantial 
improvement  on  the  precision  with  each  additional  grid 
drilled  at  200  metre  spacing  and  less# 

The  curves  plotted  in  Figure  50  are  applicable  only  to 
the  study  area  of  three  sections#  They  are  specific  to  the 
ratio  of  population  variance  to  the  population  mean  and  the 
relationship  of  the  range  of  variogram  model(s)  to  the  grid 
sample  spacings# 
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Precision  (  % )  of  average  net  coal  thickness  estimates  in  the 
small  scale  area*  Precision  has  been  calculated  at  drill 
spacings  for  three  variance  models* 


Variance  Model - 


Drill 

spacing 

Numbe  r 

of 

samples 

Statistical 

Independence 

Nugget 

Nested 

50 

3201 

.79 

.53 

.19 

100 

833 

1.54 

1.11 

.53 

200 

225 

2.96 

2.  19 

1.48 

300 

106 

4.30 

3.25 

2.70 

400 

65 

5.50 

4.24 

3.88 

600 

31 

8.01 

6.44 

6.23 

800 

21 

9.62 

8.07 

7.94 

1200 

8 

16.67 

15.  12 

15.03 

1600 

3 

28.87 

27.42 

27.34 
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4  8  10  20  40  80100  200  400  1000  2000 

number  of  observations 

Figure  50  -  Relative  precision  of  the  average  net  coal  thickness  os  a  function  of  the 
number  of  observations  in  the  small  scale  area. 
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The  equation  for  the  relative  precision  of  the 
arithmetic  mean  of  net  coal  thickness  observations  (m)  taken 
on  a  square  grid  with  spacing  d  over  an  area  A  isl 

P  =  200*S(  d  )*fT~  (5.5) 

m  \  N 

where  S(  d  )  is  the  variance  of  the  error  made  in  extending 
the  value  of  a  central  observation  to  the  enclosing  square 
cell  and  N  is  the  number  of  observations. 

Ge  nerally  f 

N  =  _A 
d2 


so 


N 


=  _42 
A 


(  5.6  ) 


Let  a  be  the  range  of  a  spherical  net  coal  thickness 
variogram  with  sill  C.  The  following  relations  holds 


S(d)  = 


d  ^  a 


(  5.7  ) 


S(  d  )  =  yj  C  d  >  a 

Substituting  (5.6)  and  (5.7)  into  (5.5) 

d^a 


P  -  200*  \Jc~  *  1  *  fd^" 
m  A  1  a 

P  =  200*  y[c~  *  d 

m  yJT 


d>  a 


(  5.8  ) 


(5.9  ) 


Three  terms  affect  the  precision  of  an  estimate.  The  first 
is  the  ratio  of  the  standard  deviation  to  the  mean.  This 
term  is  characteristic  of  the  coal  zone  and  varies  from  0.2 
to  0.6  for  examples  considered  in  the  thesis.  The  second 
term  is  the  inverse  of  the  square  root  of  the  area 


estimated.  If  the  area  studied  is  expanded  by  four  times  (at 


the  same  grid  spacing  for  observations)  then  the  uncertainty 


. 
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is  effectively  halved*  Finally  the  precision  of  the  estimate 
is  a  function  of  the  square  root  of  the  grid  spacing  cubed* 
Equation  (5*8)  may  be  rewritten  in  three  ways  to  illustrate 
these  relationships* 


P  = 

ci  \ET 

Ci  =  constant 

m 

or 

InP 

=  C2-0.5lnA 

C 2  =  constant 

or 

InP 

=  C3+1.5lnd 

C3  =  constant 

Where  a  coal  basin  has  been  explored  on  a  square  grid* 
each  coal  zone  will  be  proved  to  a  greater  or  lesser  degree 
of  precision  because  each  zone  will  have  a  characteristic 
ratio  of  standard  deviation  to  the  mean  and  a  unique  areal 
distribution*  One  uniform  grid  spacing  does  not  provide 
consistent  confidence  limits  for  all  seams,  underlining  the 
fallibility  of  arbitrary  classification  schemes  based  upon 
observation  spacing*  A  quoted  precision  is  also  dependent 
upon  the  area  estimated  and  cannot  be  applied  to  a  smaller 
area  because,  for  the  same  sample  grid,  the  precision  of  a 
smaller  area  is  higher*  Finally,  it  must  be  noted  that  these 
conclusions  apply  when  the  grid  spacing  is  less  than  or 
equal  to  the  range  of  the  major  variogram  structure  which  is 
assumed  to  have  a  spherical  form*  If  observations  are  spaced 
farther  apart  than  the  range  of  the  variogram  structure,  the 
precision  does  not  depend  upon  the  distance  and  area  but 
only  upon  the  number  of  observations  in  the  area  (formula 


5.9  ). 
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5 •  4  Fr.gsj.si9n  al£  ffalums  Estlaales 

Analytical  expressions  described  by  David  (1970)  were 
used  to  calculate  the  precision  of  the  estimated  coal 
volume*  The  relative  variance  of  the  volume  estimate  is  the 
sum  of  the  relative  variance  of  the  average  net  coal 
thlcknessf  the  relative  variance  of  the  area  estimate  and  a 
supplementary  term*  The  latter  term  is  the  product  of  the 
relative  variance  of  the  area  estimate  and  the  relative 
variance  of  net  coal  thickness  observations  and  compensates 
for  the  effect  the  uncertainty  in  the  area  has  on  the 
average  thickness  estimate*  The  precision  of  the  volume 
estimates  are  tabulated  in  Table  11*  The  precision  of  the 
average  net  coal  thickness  is  based  on  a  variogram  model 
with  a  range  of  200  metres  (  the  nested  structures  model  of 
Table  10)*  Table  11  also  shows  the  percentage  contribution 
of  the  component  variable  in  the  volume  equation  to  the 
relative  variance  of  the  volume  estimate* 

The  precision  of  the  volume  estimate  was  plotted 
against  the  number  of  observations  on  full  logarithmic  axes 
(Figure  51)*  The  area  and  average  thickness  precision  curves 
were  plotted  for  comparison*  The  volume  curve  parallels  the 
area  curve  very  closely  since  the  uncertainty  in  the  area 
estimate  comprises  over  70%  of  the  variance  for  any  of  the 
volume  estimates  (Table  6)*  The  supplementary  term  is  not 
very  significant  as  a  component  of  the  variance* 
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lafals  H 


Precision  of  the  volume  estimate  and  coaponents  of  the 
estimation  variance  at  three  drill  spacings* 


drill  spacing  (m) 

800 

400 

200 

number  of  observations 

21 

65 

225 

—  precision  - 

area 

average  thickness 

volume 

—  component  of  variance 

area 

average  thickness 
supplementary  term 


88 

71 

70 

9 

25 

27 

3 

4 

3 

25  6.5  2.4 

7.94  3.88  1.48 


26.7  7.69  2.86 


volume 


100% 


100% 


100% 
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Figure  51  -  Relative  precision  of  the  estimate  of  the 
volume  of  coal  as  a  function  of  the  number  of 


observations  in  the  small  scale  area . 
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CHAPTER  6 


Summary  and  Conclusions 


6.1  £flflclusjl<?fis  <211  the  Spatial  Cqntlnul  ty  <2l  Gaal  Zaaes 


The  continuity  of  the  thickness  and  the  structural 
surface  of  some  coal  zones  was  studied  using  the  technique 
of  structural  analysis* 


Net  coal  thickness  appears  to  be  wide— sense  stationary 
and  isotropic*  resulting  in  experimental  variograms  which 
rise  from  the  origin  to  a  range  of  influence  at  which  point 
the  variogram  is  bounded  forming  a  sill*  These  experimental 
variograms  are  commonly  composed  of  two  nested  structures:  a 
short  range  structure  and  a  long  range  structure*  The 
absolute  height  of  the  sill  and  the  relative  height  of  the 
nested  structures  varies  from  zone  to  zone* 


The  ranges  of  some  of  the  net  coal  thickness  models 
were:  200  m*  800  to  1000  m  and  4000  to  6000  m*  No  attempt 
was  made  to  relate  the  range  of  influence  to  a  geological 
model  of  coal  deposition  which  sight  explain  the 
observations  because  the  true  variability  of  the 
depositional  process  could  be  masked  by  the  size  of  the 
area*  the  location  of  the  samplesf  seam  correlation  errors* 
the  potential  mining  zone  selection  and  the  method  of 
computing  experimental  variograms* 


The  elevation  of  the  bedrock  and  the  elevation  of  the 
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coal  zone  were  non— s tat ionary  for  distances  greater  than 
1000  m  and  could  be  modelled  with  a  first-order  local  drift 
component  and  a  bounded  variogram  with  a  range  of  about  1000 
m •  The  sill  of  the  variograms  varied  for  each  individual 
coal  zone  but  ranged  from  400  ft2  to  600  ft-2*  The 
concentration  of  observations  near  the  subcrop  edge  of  the 
coal  zone  resulted  in  computed  variograms  which  represent  a 
more  structurally  disturbed  area  of  the  zone  than  the  bulk 
of  the  area  where  the  zone  is  undisturbed  by  glacial  action* 

Evidence  from  the  Estevan  coalfield  suggests  that  coal 
zones  may  be  disturbed  structurally  within  the  25  foot 
interval  below  the  bedrock  surface  and  that  the  average 
width  of  these  structures  is  less  than  800  metres  (1/2 
mile)*  The  implication  for  resource  analysis  is  that  in  the 
absence  of  drilling  at  400  m  centres  of  less*  the  location 
of  subcrop  is  uncertain  and  that  boreholes  which  intersect 
the  area  less  than  25  feet  below  bedrock  have  a  limited  area 
of  influence  for  the  purpose  of  estimating  the  area 
underlain  by  coal*  Results  from  the  small  scale  area  confirm 
that  zone  structure  may  be  locally  complex  and  that  the 
location  of  subcrop  is  not  a  simple  contour  especially  where 
an  unconformity  cuts  the  underlying  beds  at  a  low  angle* 

6.2  Conclusions  on  S_LalX&±A£_fiA  Mslfrodg 

Statistical  aethodsi  in  Combination  with  geological 
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geologic  assurance  concerning  the  estimated  quantity  of 
potential  resource  in  some  coal  zones* 

Methods  of  estimation  based  on  the  more  appealing  model 
of  spatial  correlation  (particularly  those  derived  from  the 
Theory  of  Regionalized  Variables)  resulted  in  narrower 
confidence  limits  than  methods  that  assume  statistical 
independence  of  observations*  The  former  techniques  were 
used  to  estimate  the  average  thickness  of  net  coal  in  a  zone 
over  an  area  where  sufficient  data  are  available  to  obtain  a 
stable  estimate  of  the  variogram*  Kriging  proved 
particularly  valuable  as  an  unbiased  estimator  of  averages 
over  irregularly  sampled  areas*  The  precision  of  the 
arithmetic  mean  of  observations  taken  systematically  on  a 
square  grid  was  obtained  easily  using  auxiliary  functions 
once  a  variogram  model  was  identified*  These  techniques  may 
be  of  limited  use  in  regional  coal  resource  evaluations 
because  commonly  the  data  are  not  sufficient  to  compute  a 
stable  estimate  of  the  variogram* 

The  relative  variance  of  the  areal  proportion  of  coal 
was  calculated  using  an  approximation  formula  where  samples 
were  available  on  a  square  grid  and  the  grid  cell  size  was 
smaller  than  the  average  diameter  of  eroded  areas*  The 
latter  assumption  proved  difficult  to  verify  and  limits  use 
of  the  method  to  areas  where  such  an  assumption  can  be  made 
a.  priori*  Prediction  of  the  location  of  coal  zone  subcrop  at 


the  bedrock  unconformity  using  kriging  to  estimate  surfaces 
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and  grid  manipulation  was  not  a  succes 
statistical  estimation  methods  be 
disturbance  of  the  seams  localized  with 
bedrock  surface.  The  size  of  such  feat 
the  interpolation  grid  spacing  and  smal 
spacing  over  most  of  the  subcrop  ar 
subcrop  location  using  these  methods  r 
to  be  located  closer  than  400  m  from  ea 
interpolation  grid  spacing  is  200  m  or 


sful  application  of 
cause  of  structural 
in  25  feet  of  the 
ures  was  smaller  than 
ler  than  the  sample 
ea.  Estimation  of  the 
equires  observations 
ch  other  and  that  the 
less. 


Methods  of  subjective  probability 

quantification  of  the  uncertainty  in  estimates  of 
underlain  by  a  coal  zone  where  methods  based  upon 
probability  could  not  be  used  because  of  irregular 
patterns  or  insufficient  data. 


permitted 
the  area 
ob J ec  t ive 
sampling 


6.3  Sana I a  im  geologic 
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Coal  resource  quantity  refers  here  to  an  estimate  of 
the  volume  of  coal  in  a  zone  which  is  of  sufficient 
thickness  to  be  potentially  mineable.  Areas  of  low 
overburden  thickness  or  greater  than  average  net  coal 
thickness  were  not  selected  nor  were  areas  removed  where  the 
zone  had  been  mined  or  where  mining  is  prevented  by 
overlying  cultural  or  surface  features. 
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resource  whereas  the 
classified  as  an  ind 
required  20 %  precision 
estimate  after  the  firs 
area  was  classified  as 
from  subsequent  years  w 


Boundary  zone  estimate  would  be 


icated  resource,  on  the  basis  of  a 
in  the  estimate*  The  coal  quantity 
t  year  of  drilling  in  the  small  scale 
an  Indicated  resource  while  estimates 


ere  termed  measured  resources. 


Over  two— thirds  of  the  estimation  variance  of  the 
volume  of  coal  can  be  attributed  to  uncertainty  in  the  area 
estimate*  In  the  estimates  at  all  grid  spacings  for  the 
small  scale  area  and  those  for  the  Estevan  area,  the  area 
component  of  uncertainty  was  at  least  twice  that  of  the 
average  thickness  component*  That  these  results  agree, 
arrived  at  independently  by  methods  based  on  subjective  and 
objective  probabili ty,  lends  credence  to  the  validity  of 
both  approaches* 


The  uncertainty  in  the  area  is  composed  of  the 
uncertainty  in  the  location  of  the  subcrop  perimeter  and  the 
uncertainty  in  the  location  of  lithologic  boundaries  (both 
external  and  internal)*  Although  coal— bearing  strata  of  the 
Interior  Plains  are  relatively  undisturbed  in  comparison  to 
mountain  areas,  structural  features  of  small  dimension  in 
combination  with  a  very  low  angle  of  intersection  between 
the  zone  and  the  bedrock  unconformity  cause  considerable 
uncertainty  in  the  location  of  subcrop*  Variations  in  the 
bedrock  topography  form  a  secondary  source  of  uncertainty* 
The  understanding  of  facies  relationships  in  fluvial  systems 
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has  not  yet  reached  the  state  where  the  location  and  size  of 
the  meander  belts  of  an  ancient  river  for  example)  can  be 
predicted  with  precision  given  the  subsurface  data  available 
in  regional  resource  evaluation  programs* 


Studies  have  shown  that  a  regular  square  grid  is  the 


m 


ost  efficient  sampling  pattern  for  classification  of  facies 


and  estimation  of  the  area  underlain  by  coal*  The  grid 
spacing  would  need  to  be  400  m  or  less  to  map  precisely 
lithologic  boundaries  and  the  subcrop  perimeter  of  the  coal 
zone*  Should  a  larger  grid  spacing  of  800  m  be  chosen,  the 
boreholes  must  penetrate  at  least  to  a  regional  marker  below 
the  interval  of  disturbed  bedrock  to  facilitate  zone 
correlation  and  estimation  of  regional  trend* 

The  estimate  of  the  average  thickness  of  net  coal 
contributes  the  smaller  component  to  the  uncertainty  of  the 
volume  estimate*  Geostatlst ical  methods  provide  lower 
variance  estimates  than  methods  based  on  statistical 
independence  but  the  difference  is  most  significant  when 
closely-spaced  samples  are  used  to  form  the  estimate*  The 
precision  of  an  estimate  of  average  thickness  of  net  coal 
was  seen  to  be  a  function  of  three  factors:  the  natural 


var 


iation  in  thickness  of  the  zone,  the  size  of  the  area 


over  which  thickness  was  averaged  and  the  spacing  of 
observations  made  on  a  square  grid*  For  this  reason, 
estimates  of  different  sized  areas  from  the  same  size 
observation  grid  results  in  different  relative  precision 
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values*  a  major  shortcoming  of  resource  classification 
methods  based  only  upon  distances  between  observations* 

6*4  g.gcpnHBefld&ligjas  i9£  Fvr.ttosr  Slady 

The  behavior  of  the  variogram  of  net  coal  thickness 
should  be  investigated  for  observations  located  very  close 
to  one  another*  as  this  would  contribute  to  the  modelling  of 
coal  thickness  variograms  for  data  points  from  the 
subsurface*  The  observations  could  easily  be  made  in  outcrop 
or  as  mining  progressed  along  a  box  cut  in  a  prairie  strip 
mine  • 

The  development  of  the  Theory  of  Regionalized  Variables 
was  sparked  by  problems  encountered  in  estimating  the 
average  gold  content  of  a  volume  of  rock  from  samples  of 
differing  volumes  at  varied  locations*  The  vertical  and 
lateral  variability  of  coal  quality  parameters  would  be  an 
excellent  study  topic  where  samples  are  available  from 


continuous  outcrop  or  mine  exposures* 
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